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This document contains a description of each of the software modules 
of the Image DAta Processing System (IDAPS) developed under NASA contract 
NAS8-25471 during the period September 15, 1972 to September 15, 1973. Much of 
this software is derived from an earlier effort which was reported in SDC 
TM-(L) -HU-033 /006/00 "ATM Experiment S-056 Image Data Processing System Software 
Development - Volume II". The changes which this document reflects are the 
result of additions to the application software of the system and an upgrade of 
the IBM 7094 Mod I computer to a 1301 disk storage configuration. The addition 
of the disk to the system eliminates most of the manual steps associated with 
storing intermediary data files on magnetic tape, thus greatly improving ease of 
use of the system. 

The purpose of this document is to supply necessary information 
about IDAPS software to the computer programmer who desires to make changes 
In the software system or who desires to use portions of the software outside 
of the IDAPS system. Generally, each software module is documented in the 
following manner: 

° Module Name 
® Purpose 
° Usage 

° Common Block(s) Description 
° Method (algorithm of subroutine) 

° Flow Diagram (if needed) 

° Subroutines Called 
® Storage Requirements 
° Required by: (IDAPS Operator) 
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A general description of the common blocks used by the IDAPS system 
is included as Appendix A to this document. Appendix B is an explanation of the 
diagnostic features of the system, and Appendix C describes in detail the proce- 
dure for adding new IDAPS operators to the system. 

For the convenience of the user, subroutines herein documented are 
arranged in alphabetical order. 
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Purpose : 


AREAC (Area Calculations) 

AREAC calculates the area on a sphere represented by 
pixels of a specific value. Aspect ratio is accounted 
for as a function of the distance to each pixel from the 
center of the sphere. Also calculated are the centroid 
of the area, the center of radiation intensity over the 
area, and the center of total released energy over the 

area. 


Usage: 


CALL AREAC 


Common block COM1 parameters: 


Common 

name 

Local 

name 

IN (5) 

JIN - 

SN (5) 

KIN - 

L 

NL - 

C 

NC “ 

Kl 

- r .QASE 

K2 

LC 1 - 

K3 

KC J 

K4 

LSF1 

K5 

KSFj 

K6 

NTAB 


input file containing subsetted data (output from 
SUBSET operator) 

Original input image (used with MOMENT operator only) 

Number of lines 
Number of columns 

-Pixel value for which the area is calculated 

Line and column location of the center of the sphere 
with respect to the upper left corner of the original! 

input image 

. Location of the upper left corner of the data 
frame being processed with respect to the upper 
left corner of the original input image 
- indicator of the type of aspect correction desired: 
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NTAB = 0 - No aspect corrections are applied 
NTAB = 1 “ Corrections based on table supplied by 
user 

NTAB = Other - Internal aspect ratio table used 
I0RB - Indicator of sphere 
IORB = 1 - Earth 
I ORB = Other - Sun 
IOP - Operator indicator 

IOB = 1 - AREA calculated only 
IOB = 2 - MOMENT package 

SPOT - Size in micrometers of scanning aperture used 


The centroid of a thin plate on the surface of the hemi- 
sphere is located in rectangular dimensions by: 

- EcAs - _ EgAs - _ Eft As 

X s ~ EAs Y s EAs Z s EAs 

The values of c and % can be obtained by direct methods 
and the values of g can be computed as: 



As Is the elemental area corresponding to one picture 
element and located by l , c, and g. 


For calculating the center of radiation intensity 
instead of the area, the above equations must be 
modified by replacing the elemental area s with an 
intensity factor. Thus: 



Eci 

Ei 


Y. = 

i 


Egi 

Ei 


Z 


■ _ Efcj 
i Ei 
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Likewise, for the center of total released energy, both 
the rate of energy release (intensity) and the area over 
which it is released must be taken into account so that: 


y _ ZciAs — _ IgiAs — Z&i&s 

e ZiAs Y e“ ZiAs e ZAs 

Prom the observer's point of view, the centroid of a 
feature on the sun may more easily be located by spherical 
coordinates of latitude and longitude (0 and 0) . 


0 = cos _1 (Z/R) and 6 = tan'" 1 (Y/X) 

Converting the formulas for the centroid, center of intensity, 
and center of energy; 


77 As, 


- -iip.i 

V cos ‘Hi’ 


— -l,ZAiAs, 

‘SHE’ 


— , -l,ZgAs v 

e s = tan W 


e = tan'hlai) 

i L ci 




Subroutines: DISKR, DATAIN, CONVRT 


System library routineLS,- ARCOS, SQRT, ATAN 


Storage : 


12227 (8 


words 


Required by: 


AREA, MOMENT 
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ARITH (Arithmetic Functions) 

Purpose: ARITH performs the following point by point operations 

on two frames of data: A+B, A-B, A*B, A/B. Also A+C, A*C r 
C/A and A**C, where C is a constant, may be performed. 

Usage: CALL ARITH 


Common block COMl parameters : 


Common 

name 

Local 

name 

IN (5) 

11, 12 - Input file names for A and B 

ON 

10 - Output file name 

L 

NL - Number of lines 

C 

NC - Number of columns 

K1 

KF - Operation flag 


KF = 1 - A+B 


KF = 2 - A-B 


KF = 3 - A*B 


KF = 4 - A/B 


KF = 5 - A+C 


KF = 6 - A *C 


KF = 7 - C/A 


KF - 8 - A**C 

A1 

CN - Constant to be used if KF is greater than 4 


Subroutines: DISKR, DISKW 

Storage: 14243^ words 


Required by: MATH 
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Purpose : 
Usages 


Storage : 


ARROW 


ARROW creates the image of an arrow in a 9x9 matrix. 

CALL ARROW (MAT, INDIC) 

MAT - The given 9x9 matrix 
INDIC - Arrow direction indicator 
INDIC = 1 - SW to NE 

INDIC = 2 - SE to NW 

INDIC = 3 - NW to SE 

INDIC = 4 - NE to SW 

511 words 


Required by : FRAME 
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ASSIGN (Assign Name) 

ASSIGN designates a unique name for every input, scratch, 
and output disk file for an operator. 

CALL ASSIGN (I , IN, SN ,0N ,NF ,LOF ,CARD) 

I - Operator index into array NF 
IN - Input file name array 
SN - Scratch file name array 
ON - Output file name 

NF - Operator parameter table containing the first three 

characters of each operator ' s name , and the number of 
input, scratch, and output files to be assigned for 
that operator 

LOF - Last opened output file 
CARD - Operator card image 

The operator card image is scanned for the presence of multiple 
or non-standard input file names enclosed in parentheses. If 
any are found, they are stored in array IN and replaced by blanks 
in the card image array. Otherwise, the file found in LOF is 
stored in IN(1). Standard output file names are established by 
the following convention: 

The first three characters of the operator's name are 
preceded by an integer constant which represents the 
numerical order of the file's use. For example, the first 
time operator CONVOLVE is encountered, the output file is 
named ICON; the second time CONVOLVE is encountered, the 
output file is named 2C0N, etc. 

Scratch files are named in a similar manner with the addition 
of the letter S before the operator name. For example, operator 
TRANSPOSE may have scratch files 1STRA, 2STRA, 3STRA, etc. The 
integer digit always begins at one for each operator's scratch 
files because the scratch files are deleted immediately after use. 
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i 9 


Subroutines : UNPACK , ERROR 

System library routines: CTOBXN, FLD,MOVST 


Storage : 


1305 words 
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Purpose: 


Usage: 


Method : 


Restrictions : 


BESJ (Bessel Function Generator) 

BESJ generates all Bessel functions of the first kind having 
integral order and real argument subject to the restrictions 
below. 

CALL BESJ (X, N, BJ, D, IER) 

X - Argument of the J Bessel function desired 
N - Order of the J Bessel function 
BJ - Resultant Bessel function 
D - Required accuracy 
IER - Error code 

IER = 0 - No error 

IER = 1 - N is negative 

IER = 2 - X is negative or zero 

IER = 3 - Required accuracy not obtained 

IER = 4 - Range of N compared to X not correct 

The Bessel functions of the first kind of integral order are 
computed by recurrence relation techniques described by 
H. Goldstein and R. M. Thaler, Recurrence Techniques For The 
Calculation Of Bessel Functions, MTAC, V. 13, pp. 102-108 
and A. Stegun and M. Abramowitz, Generation of Bessel Functions 
On High Speed Computers , MTAC, V. 11, 1957, pp. 255-257. 

2/3 

N < 20. + 10. x - x 
for x ^ 15. 

N < 90. + x/2. 

for x > 15. 


Storage: 


456 words 
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CENTR (Center) 

Purpose: CENTR multiplies each pixel of an input image by (-l)**(i+j) 

where i is the line position of the point and j is its column 
position. 

Usage: CALL CENTR 

Common block COMl parameters : 


Common 

name 

Local 

name 

IN (5) 

IN - Input file name 

ON 

10 - Output file name 

L 

NL - Number of lines 

C 

NC - Number of columns 

Kl 

KF - Data flag 


KF = 1 - Integer 


KF = Other - Floating point 


Restrictions: The number of columns in the image must be less than or equal 

to 1024. 

Subroutines: DISKR, DISKW 


Storage : 


2206 /Q words 


Required by: 


CENTER 
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CKFMT (Check Format) 

Purpose: CKFMT compares the format of an input file for an operator 

with the format required by that operator, and converts the 
input format if necessary. 

Usage: CALL CKFMT (KIF) 

KIF - Input file format desired 
KIF = o - Floating point 
KIF = 1 - Integer 

Common block C0M1 parameters: (See Appendix A) 

Method: The format of each input file is determined from the disk 

directory. If format conversion is necessary, a temporary input 
file is opened. The name of the original input file is saved in 
array PN, while the temporary file name replaces it in array IN. 

The converted data is stored in the temporary file and the original 
input file is closed. 

Subroutines: DISKO, DISKC, IDENT, DISKR, DISKW 

System library routine: SSWTCH 

Storage: 12246 words 

lo 
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CONV (Convolution Driver) 

CONV has five entry points which perform the following two- 
dimensional Fourier transform operations on an input image: 
CONV Convolution 
DCONV - Deconvolution 
CRLATE - Correlation 
XFORM “ FFT 
IXFORM - IFFT 

CALL COW (IP,MSIZE, IN ,NTSIZE,KFLAG, II) 

CALL DCONV (IP ,MSIZE, IN ,NTSIZE, KFLAG, II) 

CALL CKLATE (IP ,MSIZE, IN ,NTSIZE , KFLAG , II) 

CALL XFORM (IP,MSIZE,NTSIZE,NOR,II) 

CALL IXFORM (IP,MSIZE,NTSIZE,II) 

IP - Point spread function file name 

MSIZE - Size of point spread function 
IN - Input image file name 
NTSIZE - Size of output file 
KFLAG - Point spread function control flag 

KFLAG = 1 - Point spread function is being input 
in frequency domain 

KFLAG = Other ~ Point spread function is being 
input in spatial domain 
II - Scratch disk file name 

Common block C0M1 parameters: 


Purpose: 


Usage: 


Common 

name 

Local 

name 

CARD (13) 

CARDl (6) , RMAG , PHA, CARD9{5) - Input card image 

IN (5) 

IA - Input disk file names 

SN (5) 

IW, IWW, SN3 , IWWW, SN5 - Scratch file names 

ON 

10 - Output file name 

L 

NSIZE - Size of input file 
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Notes : 


Restrictions: 
Subroutines : 


Storage : 
Required by: 


If the point spread function is input in the spatial domain 
(KFLAG / 1) , an extra call is made to FRX2V to transform and 
normalize the PSF. 

The point spread function and input image arrays must be square. 
FRX2V 

3456 ^ words 

FFT, IFFT, CORRELATE, CONVOLVE, DECONVOLVE 
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Purpose : 


Usage : 


CONVRT (Convert to Degrees) 

CONVRT converts a floating point number which represents an 
angle to integer degrees, minutes and seconds. 

CALL CONVRT { HDG , 1DEG , MIN , I SEC ) 

HDG - Number to be converted 
IDEG - Degrees 
MIN - Minutes 
I SEC - Seconds 

110. _ words 


Storage : 
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Purpose : 


Usage: 


Method : 


Subroutines : 
System library 


DATAIN (Data Input) 


DATIN reads free field, tabular data from cards and stores it 
in a specified array. 

CALL DATAIN (X,N,K) 

X - Array into which data will be stored 
N - Number of data values to be input 
K - Data format 

K = 0 - Integer 
K = 1 - Floating point 

Data is input in the following general format: 

A = 9. ,8. 1,6. 3, 4. 2, 2. 7, 1.6, .8, .09, 12*0.$ 

DATAIN looks for the equal sign and isolates each number 
separated by commas until a $ is encountered at which point the 
input scan is terminated. 

Multiple values may be input by using an asterisk, as shown in 
the example above where 12*0. indicates that 12 zeros are to be 
input following the value .09. 

Integer values are input without decimal points. 

Continuation cards are required if data cannot be contained in 
columns 1 - 78 of the first card. Continuation cards must start 
in column 2 with column 1 blank. 

A check is made for the correct number of input values , as well 
as the correct format of the data. 

UNPACK, PACK, INSERT 

routines; CTOBIN, MOVST, SSWTCH 


Storage : 


1021 words 
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Required by: 


ALTER, DEPENDENT ALTER, PSF GENERATOR, FEATHER 
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Purpose : 


Usage : 


Method: 


DISKO (Disk Open) 

DXSKO has three entry points : 

DISKO - Opens a disk file and initializes the file directory 
DISKC - Closes a disk file 
DISKD - Deletes a disk file 

CALL DISKO (ID , NL, NC , KD , KU , JU) 

CALL DISKC (ID) 

CALL DISKD (ID) 

ID - File identification 
NL - Number of lines 
NC - Number of columns 
KD - Data type 

KD = 0 - Floating point 
KD = 1 - Integer (6-bit) 

KU - Disk unit number to be assigned at time file initially 
opened (0 or 1) 

JU - Disk unit number that was assigned to a previously opened 
file 

A disk file may be in one of the following states: 

Open and active - The file has been opened but 

not closed or deleted: it 

can be written or read. 

Open and inactive - The file has been opened and 

closed but not deleted: it 

must be reopened before written 
or read. 

Deleted - The file has been opened, closed, 

and deleted; it cannot be 
written or read* 

When a file is initially opened, the following entries for 
as many as 100 files are made in the file directory: 
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IDT 

NLT 

NOT 

KST 

NLPT 

NT 

KDT 


KP 


KUT 

NAF 

NOB 

LTU 

ITKADR 


File identification 

Number of lines in the input image 

Number of columns per line 

Starting track number 

Number of lines per track 

Number of tracks 

Data type 

KDT = 0 - Floating point 
KDT = 1 - 6-bit integer packed 
Pointer to buffer 

KP = 0 - Buffer closed 
KP - 1-7 - Buffers 1-7 
Unit specification (0 or 1) 

Number of assigned files 
Number of open buffers 

Last track used - LTU(l) = Unit 0, LTU (2) = Unit 1 
Track address of data currently in buffer (KP) 


In addition to the entries made in the file directory, when a file 
is initially opened a core buffer may be set aside for it. 

A buffer is required if the number of columns per line requires 
less than 233 words of storage. Thus, an integer file will use 
a buffer if a line contains less than 1393 columns since each pixel 
is 6 bits and 6 pixels are packed into a 36-bit word. The use of 
buffers reduces the number of accesses required when reading or 
writing disk, since a full track of data is read into core or 
written from core for each read or write operation. For example, 
to read or write a 128 by 128 integer image requires only 7 disk 
accesses. The following table portrays the layout of various sized 
images when stored on the disk. 
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Number of lines 
per image 

Number of columns 
per line 

Number of lines 
per track 
(integer format) 

Number of tracks 
per image 



1024 

1024 

2 

512 

512 

512 

5 

103 

256 

256 

10 

26 

128 

128 

21 

7 

64 

64 

43 

2 


Buffers are used only during an operator's execution. Each active 
file is either closed or deleted following execution of the operator. 
Input file(s) - Closed 
Scratch file(s) - Deleted 
Output file - Closed 

If an output file is to be used as the input file for a succeeding 
operation, it must be reopened. 

When a file is closed, the assigned buffer is released but the 
entries in the file directory for the closed file and the file it- 
self on disk are not changed. When a file is deleted, the buffer 
is released, the entries for the deleted file in the file directory 
are deleted, and the file is purged from the disk. 

Subroutines : IDENT 

System library routines : SSWTCH , EXIT 

Storage: 3134. words 
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Purpose : 

Usage: 

Method : 

Subroutines : 
System library 


DISKR (Disk Read) 


DISKR contains two entry points: 

DISKR - Reads a section of an image from disk 
DISKW - Writes a section of an image to disk 

CALL DISKR (ID ,LN ,LL) 

CALL DISKW (ID ,LN ,LL) 

ID - File identification 

LN - Line index. If negative, write backwards (i.e. , bottom 
to top of image) 

LL - Line to be read or written 

Pixels are stored on the disk in packed format (6 characters/ 
word) for integer files and unpacked format (1 character/word) 
for floating point files. Lines are buffered when the 
column size is less than 1393 (integer) or 233 (floating point) . 


UPL6, PL6, IDENT 

routines: SSWTCH, DTAKE, DCHEK, DPUT, EXIT 


Storage: 


11310 words 
(8 
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DRIVTP (TRNPS Driver) 

DRIVTP calculates the variable dimensions of the core buffer used 
in subroutine TRNPS , and calls TRNPS with the correct parameters. 

CALL DRIVTP (REVOLV) 

REVOLV - Rotation option indicator 
REVOLV = 1 - TRANSPOSE 
REVOLV = 2 - +90 DEGREE ROTATE 
REVOLV = 3 - -90 DEGREE ROTATE 
REVOLV = 4 - MINOR TRANSPOSE 
REVOLV = 5 - VERTICAL ROTATE 
REVOLV = 6 - HORIZONTAL ROTATE 
REVOLV - 1 - 180 DEGREE ROTATE 


Common block COM1 parameters : 


Common 

name 

Local 

name 

L 

NL - Number of lines 

C 

NC - Number of columns 

K1 

FMT - Format of the data 


FMT = 1 - Integer 


FMT = Other - Floating point 


Purpose : 


Usage: 


Method: The maximum number of columns which can be in core is determined 

by dividing a single 10240 word array by the number of lines in 
the input image. This calculated dimension, INDEX, is passed to 
TRNPS where the core buffer is dimensioned dynamically (INDEX, NL) . 

Subroutines : TRNPS 


Storage : 


24157 words 
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TRANSPOSE, +90 DEGREE ROTATE, -90 DEGREE ROTATE, MINOR TRANSPOSE, 
VERTICAL ROTATE, HORIZONTAL ROTATE, 180 DEGREE ROTATE 
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Purpose : 


Usage : 


Method : 


Restrictions : 


Subroutines : 


System library 


Storage : 


DSPLIN (Display Line) 

DSPLIN converts a string of BCD characters into a 9x171 array for 
display of alphanumeric characters on the Photowriter or DICOMED 
display unit. 

CALL DSPLIN (IN , NC , KC , KOUT) 

IN - String of characters to be converted to display format 

NC - Number of characters in string to be converted, 

excluding trailing blanks 
KC - Starting column position in output array 
KOUT - Output array (9x171 packed format) 

Each BCD character in the input string is interpreted and a 
corresponding 9x7 bit configuration is placed in the output array. 

A maximum of 102 characters can be converted for display at one 
time. 

UNPAC 

routines : FLD 

22566 words 
(8 


Required by: 


LABEL 
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Purpose: 


Usage: 


Method : 


Restrictions : 


System library 
Storage: 


DSPLN (Display One Line at a Time) 

DSPLN converts a string of BCD characters into an array for display 
of alphanumeric characters on the Photowriter or DICOMED display. 
This subroutine is used for the line-at-a-time output of labels 
required by the FRAME and MULTIPLE DISPLAY operators. 

CALL DSPLN (IN ,NC , KC ,OUT ,TP ,NUM ,LN) 

IN - String of characters to be converted for display 
in BCD format 

NC - Number of characters in string to be 
converted 

KC - Starting column position in output array 
OUT - Output array unpacked 

TP - Maximum number of columns in the output array 
NUM - Index to first character in the input array 
LN - Line number (1-9) 

Each BCD character in the input string is interpreted and a 
corresponding bit configuration is placed in the output array. 

DSPLN must be called in a loop, 9 times, with LN incremented each 
time, to effect the correct 9x7 bit character configuration. 

A maximum of 80 characters can be converted for display at one 
time. 

routines: FLD 

436 words 


Required by: 


FRAME, MULTIPLE DISPLAY 
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DSTANC (Distance) 

Purposes DSTANC calculates the distance on a sphere from one point to another, 

given the latitude and longitude of each point. 

Usages CALL DSTANC 

Common block C0M1 parameters: (See Appendix A.) 

Method : 


A point (A) is located on the surface of a hemisphere of radius R. 
Two sets of coordinates are chosen to describe the iocation of A - 
spherical coordinates 0, <{> , and R, and rectangular coordinate 
angles a , t> ; ' and y . 
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= The angle that OA makes with the 2 axis. 

0 = The angle between the projection of OA on the XY plane and 
the Y axis . 

a = The angle between OA and the X axis. 

3 = The angle between OA and the Y axis, 

y = The angle between OA and the 2 axis. 

If the sphere is viewed at a great distance by an observer situated 

on the Y axis, then all points on the hemisphere will appear as a 
perpendicular projection on the X/Z plane. All lines which are 
parallel to a radius vector in the X/Z plane will appear in true 
perspective while all others will appear foreshortened by some 
aspect ratio. 


From the viewer ' s observations , the following calculations can be 
made : 

R = 1/2 the diameter of the apparent solar disk 
<j> = cos (— ) 

. -1 , c , . -1 f c 

6 = sin ( Oo ." - J ® sin Tty 

Rsin * {/r 2 i 2 } 

-1 C 


a ~ 


cut, R 


, -1 y/c 2 + 9? 


~ sin 


* -1 ft 

* - cos - 


-1 


a - cos (sine}) sinB) cf> = y 


6 = 


-1 


-1 ,cosS, 


cos (sir4 cos0) and 0 = cos ( — : ) 

T siny 


Y = 4> 


For two vectors 0 — and 0 — the cosine of the angle between 
then is: 


cos 



cosa^ cosa^ + cosg^ cosg^ + cosy^ cosy 2 
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Translating this to spherical coordinates: 


OA 


L I -1 

0A = cos (sin^ sin<t> 2 sinO^ sin0 2 + sin^ sin<}>. 


cosQ^ cosQ^ + cos<f>^ cos<f> 2 ) or 


-1 


= cos (sinij) sin<f> (sin0. sin0 + cos9 cos6 ) + cos<{> cosifO 


since, sinQ^ sin9 2 + cos ^i cos ®2 = cos ^1 “ ® 2 } 
Then : 


OH 


/_ 0H = cos --1 (sin<j> 1 sin4> 2 (cos^ - 6^) ) + cos^ cos^) 


and the distance between OA^ and 0A 2 is R 


OA. 

/ 1 

0A„ 


Subroutines : 


CONVRT 


System library routines: CTOBIN, SIN, COS, ARCOS, SQRT 


Storage : 


514 words 
(8 


Required by: 


DISTANCE 



15 September 1973 


29 


System Development Corporation 
TM -(L) -HU-033/008/00 


DUMPF (Dump File) 

Purpose: DUMPF has two entry points: 

DUMPF - Dumps a disk file to tape 
READF - Loads a file from tape to disk 

Usage: CALL DUMPF 

CALL READF 

Common block COMl parameters: 


Common 

name 

Local 

name 

IN (5) 

IN - Input disk file names 

ON 

ON - Output disk file name 

L 

NL - Number of lines 

C 

NC - Number of columns 

Kl 

KT - Tape unit number 

K2 

KF - File number 


Subroutines: DISKR, DISKW 

System library routines : SKFBIN 

Storage: 4275 words 


Required by: 


DUMP FILE, LOAD FILE 
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DUMPFD (Dump File Directory) 

Purpose: DUMPFD saves a current copy of the disk file directory on the disk. 

Usage: CALL DUMPFD (LOF ,NF) 

LOF - Last opened output file 
NF - Operator parameter table 

System library routines: DPUT, DCHEK 


Storage : 


274 3, „ words 
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Purpose ; 


Usage: 


Storage: 


ERROR 

ERROR causes the printing of a message on the on-line printer . 
After a pause, the next card is read in from the card reader. 


CALL ERROR (CARD) 

CARD - Array of 13 words containing the packed (6 ch/wd) 
card image 


103 (8 


words 
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FFIL (Fourier Filter) 


Purpose: FFIL filters an input image in the frequency domain. 

Usage: CALL FFIL 

Common block COM1 parameters: 


Common 

name 

Local 

name 

IN (5) 

IN (1) - Image file name 


IN (2) - Filter file name 

SN (5) 

SN - Scratch files 1 through 5 

ON 

ON - Output file 

L 

NL - Number of lines 

i 

C 

NC - Number of columns 


Method: The following sequence of operations is performed: 

1) The image is transformed using subroutine FFT. 

2) The transform is converted from real/imaginary to 
magnitude/phase format. 

3) The magnitude terms of the transform are multiplied by 
corresponding terms of the filter. 

4) The resulting transform in magnitude/phase is converted 
back to real/imaginary. 

5) The transform is inverse transformed. 

1) The filter must be the same size as the image. 

2) The image and filter must be square. 


Restrictions : 
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Subroutines ; 


Storage : 


FRX2V, DISKR, DISKW 
3362 _ words 


Required by: 


FOURIER FILTER 
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Purpose: 


Usage : 


Storage : 
Required 


FLIP (Vertical Rotate) 


FLIP inverts the order of the elements in a line of data 
so that column NC becomes column 1, column (NC-1) becomes 
column 2, etc. 

CALL FLIP (NC,IBUF) 

NC - Number of columns in a line 
IBUF - Input buffer for one line 

2070 words 
vo 

by: TRANSPOSE, +90 DEGREE ROTATE, -90 DEGREE ROTATE, MINOR TRANSPOSE, 

HORIZONTAL ROTATE, VERTICAL ROTATE, 180 DEGREE ROTATE 
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FLTOFX (Floating to Fixed Point Conversion) 

Purpose: FLTOFX converts an input image from floating point to fixed point 

format . 

Usage : CALL FLTOFX 


Common block COM1 parameters : 


Common 

name 

Local 

name 

IN (5) 

IK - Input file name 

ON 

IO - Output file name 

L 

NL - Number of lines 

C 

NC - Number of columns 


Method: Each pixel within the image is rounded to the nearest whole number 

and is output in integer format. 

Subroutines: DISKR, DISKW 

Storage: 10147 words 

\8 _ - - 


Required by: 


FIX 
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FRAME (Standardized Frame Mask) 


Purpose: FRAME provides a standard frame for a processed image 

consisting of the following information: 

1) Gray scale step wedge 

2) Tic marks along the edges of the image 

3) Optional grid 

4) Optional labels, both above and below the image 

5) Histogram of gray scale distribution 

6) Border around the image 

7) Optional arrow to point at an area of interest 

8) Optional inversion of gray scale values 

Usage : CALL FRAME 

Common block C0M1 parameters: 


Common 

name 

Local 

name 

CARD (13) 

KARD (13) 

IN (5) 

ID - Input file 

L 

TNL - Number of lines 

C 

TNC - Number of columns 

Kl 

TL - Upper label 

TL = 0 or blank - No top label 

TL > 0 - Number of lines in top label 

K2 

BL - Bottom label 

BL - 0 or blank - No bottom label 

BL > 0 - Number of lines in bottom label 

K3 

G - Grid option 

G - 0 or blank - No grid 

G = 1 - Automatic grid 

G > 1 - User supplied grid increments 

K4 

OPl - Gray scale inversion option which inverts the image only 
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K5 

K6 

K7 

K8 


Method : 


Restrictions: 


Subroutines : 
Storage: 
Required by: 


0P2 - Gray scale inversion option which inverts the entire 
output frame 

AL - Arrow point line coordinate 
AC - Arrow point column coordinate 

INDIC - Arrow direction indicator. If 0, then SW to NE direction 
is used unless image borders indicate otherwise. 

INDIC = 1 - SE to NE 
INDIC = 2 - SE to NW 
INDIC = 3 - NW to SE 
INDIC = 4 - NE to SW 


The designated information is overlayed on the input image, one 
line at a time , and written onto an output tape with the gray scale 
vredge lined up vertically along the right side of the image and the 
histogram below all specified bottom labels. The completed frame 
width will be TNC + 106. 

(1) The input image may not exceed 2048 x 2048. 

(2) Only one arrow may be specified. 

(3) The labels are centered, one 80 column card image per line, 
over the entire width of the framed image. 

(4) The label may not exceed (TNC + 106) /10 characters. 

FRAMWR, HISTOW, ARROW, DISKR, DSPLN 

7022 (q words 

FRAME 
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Purpose : 


Usage: 


Subroutines : 


System library 
Storage : 


FRAMWR (Frame Write) 


FRAMWR packs a buffer and writes it on the output tape, A6, 
a designated number of times (NUM) . If the gray scale inversion 
option is indicated, the inversion is performed before the record 
is written. 


CALL FRAMWR (RECORD ,LEN , IOP , NUM , PL , PACKED) 

RECORD - The record to be written out on tape 
LEN - Length of the record 

IOP - Gray scale inversion option indicator 
IOP = 0 - No inversion 
IOP = Other - Invert 

NUM - The number of times this record is to be written 

PL - The packed line buffer 

PACKED - The length of the packed buffer 


PL6 

routines: WRITER 

4127 ^ words 


Required by: 


FRAME 
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Purpose : 


Usage : 


Notes : 


Method: 


FRXFM (Fast Fourier Transform Routine) 

FRXFM computes the Fast Fourier Transform of a complex vector 
using a Radix 4 plus 2 algorithm. 

CALL FRXFM (N2P0W,X,Y) 

N2P0W - Log 2 N, where N is the number of elements in the 
complex array 

X - Array containing the real part of the complex vector 

Y - Array containing the imaginary part of the complex 

vector 

The array need be dimensioned only in the calling routine. If a 
real array is to be transformed, set the Y array elements to zero 
and proceed as usual. 

The transform is performed in place; therefore, the input data is 
destroyed as the transform is computed. 

The routine is a modified implementation of a procedure originally 
developed by Cooley and Tukey and later improved by Gentleman and 
Sande . 

If N is the size of the vector to be transformed, then the equations 
describing discrete Fourier transformations are given by: 



matrix form as 


then equation (1) can be written in 
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2 ... ... 

In general, equation (4) requires N multiplications and divisions. 

If, however, N has M factors r , r , r such that n = H r., 

± ; m i=i 

then T can be decomposed into the product of m elementary trans- 
formations , followed by a permutation P of the result: 

T = P T T , T. (5) 

m m-1 1 


Each step T. is in turn composed of N/r. transforms of dimension r . 

3 3 3 

Since the number of arithmetic operations for a tranform of dimen- 

2 

sion r . is of the order r . , the total number of operations for 

^ m 

v in 

the transformation T is less than N I r . = mNr but N = r or 

3=1 ° 

log N-m, therefore computation requires r^ log^N arithmetic 
r 2 

operations compared with the N operations required by a standard 
transform. Also by reduction of arithmetic operations, accuracy 
is improved. 

The recursive equations needed to implement the Fast Fourier 
Transformation routine for arbitrary length records (r = 2 a for all 
integer a) were derived by G. D. Bergland and are shown below. 



15 September 1973 


41 


System Development Corporation 
TM- (L) -HU-033/008/00 


A P ( V j l'”' ,:5 p-l' k n-p-l"*‘ 'V 

r -1 
P 

= -l 

k • =0 * * 

n-p 

W N tj p-l (r i r 2*** r P -l ) + --‘ + j 0 ] k m-p (r P +l*‘* r m ) (6) 


The last entry thus calculated will yield the Fourier sum as 


X(j 


m-1 


f • ♦ • f 3 q) = 


VV 


■w 


(7) 


in such an order that the index of an X must have its binary bits 
put in reverse order to yield its index in the array 


In some cases the computation required to evaluate equation (6) 
can be reduced by regrouping the equation as 


A p { VV*’" j p-l' k m-p-l'"‘ ,k 0 ) 

fr -1 i k 

P ^p-1 m-p 

.. I Vl (j 0 ,j l j p-2' Vp'-'-'V w r 

K — U lr 

h n-p 

„ j , (k _ (r_...r ) + ...+ ^k+r) (rr .r ) 

W J p-1 n-p-1 p+2 m 1 n 0 1 2 p-l 18J 

N 


The bracketed term represents a set of r^-point Fourier transforms 
and the complex exponential weights outside the brackets reference 
each set of results to a common time origin. 


The term 


W 

r 

P 



2iri/r 
e p 


forms the basis for the complex exponential weights required in 

evaluating each r -point transform, and j . , k are the indices 
p p — x m — p 

of the transform. 
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If the , r 2 '***' r m are ec J ua - ! - to some integer, sa ¥ a > ttlen t * ie 

transformation described by the resulting recursive equations is 

referred to as a Radix a transformation. Cooley and Tukey stated 

in their original paper that the FFT algorithm is formally most 

efficient when the number of samples in a record is a power of 3 

(i e . r =3) but such an algorithm is only slightly more efficient 
j 

than a Radix 2 or a Radix 4 routine. Later it was recognized that 
the symmetries of the sine and cosine weighting functions made the 
Radix 4 algorithm more efficient than either the Radix 2 or Radix 3 
algorithm. Radix 4 algorithms, however, suffer from the severe 
limitation on acceptable record sizes. To^ overcome this difficulty 
Gentleman and Sande developed an algorithm which performs as many 
iterations of the transform as possible with a Radix 4 routine, 
and then, if required, performs the last iteration with a Radix 2 
routine. A modified version of Gentleman and Sande' s algorithm 
is implemented as the subroutine FRXFM. 

If N is the size of the array such that N = 2 N2POW , then N2P0W is 
divided by two (integer division) to compute the required number 
of transformations of dimension four (N4P0W) . 

Restrictions: The dimension (N) of the input array must be a power of two. if it 
is not, zeros -iuay be used to_buffer it up to a length which is 
a power of two. 

System library routines: COS, SIN 

Storage : 1 2 24 ( 8 word S 


Required by: 


FFT, IFFT, CORRELATE, CONVOLVE, DECONVOLVE, FOURIER FILTER 
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Purpose : 


Usage: 


FRX2V (Two-dimensional Fast Fourier Process Driver) 

FRX2V performs a two-dimensional Fourier transformation by 
calling the fast Fourier subroutine FRXFM first by rows then 
by columns . FRX2V also serves as the main processor for 
the convolution, deconvolution, correlation, transform, 
normalized transform, and inverse transform processes. 


CALL FRX2V (NSIZE, NPSIZE, IN, IW, IO, IM, IAPPLY, IP, II, 
RMAG, PHA, ITRUNC, TMAX, TMIN) 

NSIZE - Input image dimension 
NPSIZE - Output image dimension 
IN - Input image disk file name 
IW - Scratch disk file name 

10 - Output image disk file name 

IM - Forward/Inverse transformation flag 
IM = 0 - Forward transformation 
IM = 1 - inverse transformation 
IAPPLY - Process control flag 
IAPPLY = 1 - Convolve 
IAPPLY = 2 - Deconvolve 
IAPPLY = 3 - Transform and normalize 
IAPPLY = 4 - Correlate 
IAPPLY— = Other - Transform _ 

JP - - Transformed point spread function disk 
file name 

11 - Scratch disk file name 

RMAG “ Limit on magnitude amplification — See SPRED2 

PHA - Limit on phase amplification — See SPRED2 

ITRUNC - Truncation flag 

ITRUNC = 1 - No truncation 
ITRUNC - 2 - Truncate values outside of 
TMIN and TMAX 

ITRUNC = Other - Truncate values outside of 
0 . and 63 . 
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Method : 


TMAX - Maximum truncation limit 
TMIN - Minimum truncation limit 


A two-dimensional Fourier transformation is achieved by 
first performing a one-dimensional transformation of each 
of the rows of the original array and then performing a 
one-dimensional transformation of each of the columns of 
the array resulting from the first set of transformations. 
To better understand this procedure consider the continuous 
two-dimensional transform of equation (1) . 


h(u,v) - rrkcx.Y) e- 12 " ux e- I2,, V dxdy (1, 

--CO - '— as 

If an Interval T is divided into N equal parts both horizontally 
and vertically such that NAX = NAY = T, then the maximum 
frequency F = — ~ Tv' Th e incremental frequencies are then 


AU = AV = F/N so that AUAX = AVAY = — . Using the indicies j,k, 
£,tn equation (1) in the discrete case becomes 


N-l N-l 

h (£Au, mAu) = £ . £ H(jAx,kAy) e 

k=0 i*o 


-2Tii -2iri , 

j£ — — km 


N 


N 


( 2 ) 


Let 


N-l 

TR(j,m) = £ H(jAx.kAy) e 

k=0 


-2iri. 
N 


km 


(3) 


then 


N-l 


h(£Au,mAu) = £ TR(j ,ra) e 

j-o 




N 


(4) 


but a one-dimensional discrete Fourier transform over an 
Interval T, sampled N times (NAX = T and F = AU - F/N) 

is given by 

„ -2mi , 

N-l • — ^ — km 

T (MAX) = £ V(K) e 

k=0 
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Therefore, equation (3) is just N ordinary one-dimensional 
transforms, one for each row of the picture, and equation (4) 
is N ordinary one-dimensional transforms, one for each 
of the columns resulting from the N applications of 
equation (3) . 

Image buffering: If the input picture (s) and/or point spread function, 

have dimensions less than the smallest power of two 
greater than NPSIZE they are buffered with zeros equally 
on all sides up to this power. Indeed, for convolution, 
deconvolution, or correlation this is necessary in order 
to prevent image overlapping due to the cyclic nature of 
the Past Fourier Transform. In general if M is the size 
of the input picture and L is the size of the point spread 
function, the NPSIZE should be at least the smallest power 
of two greater than M + L. 

Point spread function 

transformation: If the PSF is input in the spatial domain, an extra call 

to FRX2V is made by CONV to transform and normalize the PSF 
(IAPPLY = 3) . The normalization consists of summing the 
gray levels of all elements of the PSF and dividing each 
frequency element (both real and imaginary) by this amount. 
This procedure insures than a unit amount of light will 
be dispensed for each unit present in the input picture. 

Image Centering: The transform of the PSF (real and imaginary) is also 

multiplied by (-1) 1+ ' 3 , i=l, 2,..., NPSIZE: j=l, 2,..., 
NPSIZE, for convolution and deconvolution so that the output 
picture will be centered in the frame. 

Clipping: Due to improper amplification of certain frequency terms 

due to noise, improper PSF, etc., the gray levels obtained 
from deconvolution may extend over a large dynamic range. 

The subroutine will determine the maximum and minimum 
values encountered and will print these on-line. If the 
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Data storage: 


Restrictions : 


Subroutines : 


Storage : 
Required by: 


64 levels were divided over this large range one would 
obtain an output image with little contrast. A much 
improved image will result if the output is clipped both 
above and below, with the remaining levels quantized to 
64 gray levels. Presently the subroutine clips at 0. 
and 63.? that is, all values below zero are set to 0. 
and all values above 63 are set to 63. The remaining 
range is truncated to 64 gray levels. 

The figure on the next page represents the data storage 
layout for a 128 by 128 frame of data. Since the whole 
picture will not fit into core at the same time, the frame 
is divided into four sections and each section processed 
in turn. Frames of data which are 256 by 256, and 512 
by 512 are similarly handled to minimize data transfer. 

1) The image size must be greater than or equal to 16 
by 16 and must be less than or equal to 512 by 512. 

This should be considered the maximum output image size. 

2) The output image dimension must be square and must be 
greater than or equal to NSIZE. 

3) The inverse operator assumes the input image is square. 

DISKR, DISKW, FRXFM, SPREDl , SPKED2 , SPRED3 

25, 520. words 
(o 

FFT, IFFT, CORRELATE, CONVOLVE, DECONVOLVE, FOURIER FILTER 
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*1 

Block of 


X 

.32_jw.QCis_ 

Line 1 columns 1-32 ' 

33 


Line 2 columns 1-32 

65 1 






993 


Line 32 columns 1-32 

/ 

1025 


Line 1 columns 1-32 ' 

1057 






2017 


Line 32 columns 1-32 

2049 

' 

Line 1 columns 33-64 

2081 









3041 

» 

Line 32 columns 33-64 

4 

3073 

Ua/VT'- 

\ 

Line 1 columns 33-64 1 

1 

\ 


| y\r 

1 j 

4065 

l 

L 

Line 32 columns 33-64 ] 

s 

409 7 


Line 1 columns 65-96 

4129 


etc. 


A y\y\A 



Real 


Imaginary 


Real 


Imaginary 


Real 


128 x 128 Data Storage Layout 
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Compute 

Program 

Variables 


Transform 
•v Mode? - 


forward 


Read image 
input file 
into IN 


| Inverse 

Read Image 
Transform Input 
file into IN 


Compute 

normalization 

factor 


Compute 

complex 

conjugate 


Buffer image 
size to a 
power of two 


Perform row- 
by-row 

transformatior 


Perform row- 
by-row 

transformation 


Data stored 
on scratch 
files IW, II 


Perform column- 

by-column 

transformation 


Transform 
v mode? / 


Forward 


Inverse 


Compute 
minimum & 
maximum 

r ; 

Clip Output 
if X>63, X=63 
if X< 0, X=0 


=> 

Convolution 

(complex 

multiplication) 

=1 ^ 

Deconvolution 
(complex — 


IAPPLY? 



division) 


Normalization 

(multiply 

by (-1) i+J) 


Correlation 
(complex multi- 
plication com- 
plex conjugate) 


Subroutine FRX2V Simplified Flow 


Write 

Output 

file 
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FSIZE (File Size) 

Purpose: FSIZE calculates the size of the scratch file needed by the 

two-dimensional FFT operator, according to the following 
table: 



Usage: CALL FSIZE (NPSIZE, NLINES, NCOLMS) 

NPSIZE - Output Image Size 

NLINES - Number of lines in scratch file 

NCOLMS - Number of columns in scratch file 

Storage: 112 words 

1 ° 
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Purpose : 


Usage: 


FUNGN (Function Generator) 

FUNGN produces a symmetrical function in a square array 
from an input curve. 

i 

CALL FUNGN 


V 

Common block C0M1 parameters: 


Common 

Name 

Local 

Name 

ON 

IOUT - Output file name 

L 

L - Number of lines 

C 

C - Number of columns 

Kl 

Kl ~ Input option 


Kl = 1 - Input three horizontal line segments 


Kl = Other - Input curve profile 

— 


Method: The distance from the center of the array is computed for each 

point in a square array. The distance is used to interpolate 
a value from the input curve, which is placed into the square 
array. For distance values outside the range of the input curve 
corners, the point in the square array is set to the last value 
found in the curve. 

Subroutines: DATAIN, DISKW 

System library routines: LAGRNG, COS, SQRT 

Storage: 10644 words 


Required by: 


FUNCTION GENERATOR 
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FXTOFL (Fixed to Floating Point Conversion) 

Purpose: FXTOFL converts an input image from fixed point format to 

floating point format. 

Usage: CALL FXTOFL 

Common block COMl parameters: 


Common 

Local 

name 

name 

IN (5) 

IN - Input file name 

ON 

10 - Output file name 

L 

NL - Number of lines 

C 

NC - Number of columns 


Subroutines : 


DISKR, DISKW 


Storage: 


10150 words 


Required by: 


FLOAT 
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Purpose : 


Usage: 


Subroutines: 
Storage : 


HISTOW (Histogram Write) 


HISTOW generates a histogram from an input array of gray scale 
totals for output by the FRAME operator. 

CALL HISTOW (GSD, TP, NQ, IOP, PL, PACKED) 

GSD Array of totals for the 64 gray scale values 
TP - Total number of columns in the final frame 

NQ - Buffer used for generating the gray scale indicators 
and tics 

IOP - Gray scale inversion option indicator 
PL - Packed line buffer 
PACKED - Length of PL 

FRAMWR, DSPLN 
675 , n words 


Required by: 


FRAME 
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HRDCPY (Hardcopy Printout of Images) 

Purpose: HRDCPY prints a frame of image data in 1) fixed or floating 

point format, 2) sixteen levels of gray, or 3) a map of values 
above a specified threshold* 

Usage: CALL HRDCPY 


Common block CQM1 parameters: 


Common 

name 


Local 

name 


IN ( 5 ) 

L 

C 

K1 


K2 

K3 

K4 

K5 


JIN - Input file name 
NL - Number of lines 
NC - Number of columns 
NO - Control flag 

NO - -4 - Print gray level values on-line in floating point 

NO = -3 - Print gray level values off-line in floating point 

NO — -2 - Prints gray level values on-line in integer format 

NO = -1 - Prints gray level vaues off-line in integer format 

NO = 0 to 63 - Prints blanks for values less than or equal 
to NO and X for values greater than NO 
NO = 64 - Prints 16 levels of gray using overstrike method 
ITOP - Starting line number 
JLFT - Starting column number 
LINES - Number of lines to be printed 
KLMNS - Number of columns to be printed 
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Method: 


Subroutines : 


Storage : 


When printing gray level output, HRDCPY divides each gray level 
value (integer values 0 to 63) by four, adds one and then prints 
up to nine overstrikes of selected characters to produce one 
of sixteen different gray levels. 

DISKR 

15104 words 
(o 


Required by: 


PRINT 



15 September 1973 


55 


System Development Corporation 

TM- (L) -HU-033/008/00 


ID APS 2 (Executive Control For IDAPS ) 

Purpose: IDAPS2 is the main control program for the system and performs, 

or controls the performance of, the following functions: 

1) Initialization of the system 

2) Input and interpretation of IDAPS operators and their 
associated parameters 

3) Assignment of input, output, and scratch file names 

4) Calling the necessary subroutines to execute the operator 

5) Deletion of the assigned scratch files at the end of 
execution of each operator 

6) Closing of input and output files 

Common block COMl parameters: (See Appendix A.) 

Subroutines: IN IT, SCK, OP1-OP66, OPA-OPH, XPAND, SHADE, PDGSA, MLTAVS, 

DCONV, RS1GEN, HEDCPY, PSFGEN, CONV, XFORM, IXFORM, CRLATE , 
SHAVE, FXTOFL, FLTOFX, SPINT, VARY, MAC, LABL, CENTR , MERG, 
SPLT , RITOMP, MPTORI, ARITH, SCAL, NOSFUG, INIMAG, OTIMAG, 
STKPIC , MULDSP , FRAME, REWND, DRIVTP, SEGMNT, DUMPF, READF , 
ISOPLT, ORAND, DELCLS , STAT, INSRT, RESTOR, SHRINK, AREAC , 
TFGEN , DSTANC , FUNGI'?, KONVL, SLICE, FFIL, NSA, TIME, DUMPFD, 
DISKD, RDC 


Storage : 


4572 


(8 


words 
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Purpose : 


Usage: 


Storage: 


IDENT (Identify File Index) 

IDENT is a function v/hich returns the index to the disk file 
directory for a specified disk file name, 

INDEX = IDENT (ID) 

ID - Disk file name 

INDEX - Pointer to disk file directory 

1663 words 
(8 
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INIMAG (Input Image ) 

Purpose: INIMAG has two entry points: 

INIMAG - Reads an input image in scanner format from 
tape unit A5. 

OTIMAG - Writes a file onto tape unit A6 in display format. 

Usage : CALL INIMAG 

CALL OTIMAG 


Common block COMl parameters: 


Common 

name 

Local 

name 

IN (5) 

IN - Input file name 

ON 

ON - Output file name 

L 

L - Number of lines 

C 

C - Number of columns 

Kl 

KF - File number on input tape (for entry INIMAG) 
Complement flag (for entry OTIMAG) 

KF = 1 - Complement 

KF = Other - No complement 

K2 

SL - Starting line 

K3 

SC - Starting column 


Subroutines: UPL6, DISKW, DISKR, PL6 

System library routines: SKFBIN, SKRBIN , R.EDTPR, WRITER 

Storage: 5046 (8 words 


Required by: 


INPUT IMAGE, OUTPUT IMAGE 
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INIT (Initialize Operator) 


Purpose : 


Usages 


INIT contains the following entry points: 

INIT - Initializes 
RDC - Read card 

DELCLS - Delete/Close disk files 

OPA-OP66 - Individual operator parameter and disk file 
manipulation 

Purpose 

CALL INIT ( I , LOF , NF , NO , KREP ) System Initialization 

CALL RDC ( I, LOF, NF, NO, KREP) Read and interpret operator card 

CALL DELCLS <I,INX,SNX,PNX,NF) Delete and/or close disk files 


CALL OPA (J,K) 

Initializes 

ALTER, FIX, FLOAT 

CALL OPB (J) 

ii 

DEPENDENT ALTER, FEATH: 

CALL OPC (J) 

it 

BAR CHART, PSF GENE RAT* 

CALL OPD 

ii 

R/I TO M/P, M/P TO R/I 

CALL OP1 

u 

SIZE 

CALL OP 2 

it 

EXPAND 

CALL OP 5 

it 

AVERAGE 

CALL OP 6 

ii 

DECONVOLVE 

CALL OP 7 

ii 

STOP 

CALL ..OP 9 

n 

PRINT . . 

CALL OP 11 

ti 

EXTRACT 

CALL OP 12 

n 

CONVOLVE 

CALL OP13 

ii 

FFT 

CALL OP 14 

ti 

IFFT 

CALL OP 15 

ii 

CORRELATE 

CALL OP 19 

tl 

SMOOTH 

CALL 0P20 

11 

COMMENT 

CALL 0P24 

11 

LABEL 

CALL OP25 

If 

CENTER TRANSFORM 

CALL OP26 

It 

MERGE 

CALL OP 30 

11 

MATH 
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CALL 

OP31 

Initializes 

SCALE 

CALL 

OP33 

It 

NON-SYMMETRICAL PSF 

CALL 

OP34 

II 

INPUT IMAGE 

CALL 

OP35 

11 

OUTPUT IMAGE 

CALL 

OP36 

11 

STACK 

CALL 

OP37 

11 

MULTIPLE DISPLAY 

CALL 

OP38 

ff 

FRAME 

CALL 

OPF 

11 

TRANSPOSE, +90 ROTATE, -90 ROTATE, 
MINOR TRANSPOSE 

CALL 

OPG 

11 

VERTICAL ROTATE, HORIZONTAL ROTATE 
180 ROTATE 

CALL 

OP47 

ff 

SUBSET 

CALL 

OP48 

II 

DUMP FILE 

CALL 

OP49 

11 

LOAD FILE 

CALL 

OP50 

(KQ) 

ISOGRAM 

CALL 

OPH 

II 

DELETE FILES 

CALL 

OP 53 

(NF, ID, NO) 

INCLUSIVE OR, AND 

CALL 

OP54 

It 

HISTOGRAM 

CALL 

OP55 

(KREP) 

REPEAT 

CALL 

OP56 

II 

INSERT 

CALL 

OP57 

(KQ1) 

RESTORE 

CALL 

OP58 

II 

SHRINK 

CALL 

OP59 

M 

AREA, MOMENT 

CALL 

0P61 

(A,B,E,D) 

MTF 

CALL 

OP63 

IT 

FUNCTION GENERATOR 

CALL 

OP64 

11 

FILTER 

CALL 

OP65 

It 

SLICE 

CALL 

OP66 

II 

FOURIER FILTER 


I - Operator Index 



LOF - 

* Last opened output 

file name 


NF - 

Operator parameter table 


NO - 

Number of operators 

in NF 


KREP 

- Repeat flag 



INX - 

- Input file table 



SNX - 

■ Scratch file table 



ONX - 

- Output file name 
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A,B,E ,D - Independent variable parameters for MTF operator 
J - Format desired in files 
K - Variable format 

KQ - Logical variable used with ISOGRAM operator 

JD - Array used with DELETE FILES operator 

KQ1 - Logical variable used with RESTORE operator 

Common block COM1 parameters: (See Appendix A.) 

Method: The entry for each individual operator handles reading any 

parameters into the block COMl , and the opening of all 
required input, scratch, and output disk files. 

Subroutines: ASSIGN, CKFMT , FSIZE, TSFSZ, DISKD, DISKO, IDENT, PRERUN, 

DISKC , IPARM. 

System library routines: SSWTCH, STRTMR, FLD, MOVST, CTOBIN 


Storage : 


11342 words 
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INSERT (Insert Character ) 

Purpose: INSERT is used by subroutine DATAIN to insert a character into 

a string of characters. 

Usage: CALL INSERT (CHAR, STRING, POS) 

CHAR - Character to be inserted 

STRING - String into which character is to be inserted 
POS - Position in string to insert character 

System library routines: FLD 

Storage: 100 jg words 
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INSRT (Insert Array) 

Purpose: INSRT inserts a sub-frame into a specified position in a 

larger file. 

Usage: CALL INSRT 

Common block COM1 parameters: 


Common 

Name 

Local 

Name 

IN(5) 

INI - Name of larger disk file 

IN 2 - Name of disk file to be inserted 

ON 

ON - Output disk file name 

L 

NL - Number of lines in large file 

C 

NC - Number of columns in large file 

Kl 

KL - Number of lines in sub-frame 

K2 

KC - Number of columns in sub-frame 

K3 

SL - Starting line to start insertion 

K4 

SC - Starting column to start insertion 

K5 

KF - Edge modification flag 
KF = 1 - Smooth edges 
KF = 2 - Border around sub-frame 
KF = Other - No modification 


Subroutines: DISKR, DISKW 

Storage: 21016 * words 


Required by: 


INSERT 
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IPARAM (Parameter Check) 


Purpose: IPARAM is a function which checks the format of all parameters 

on an input operator card. If the formats are correct, the 
parameters are converted to binary and stored in the common 
parameter area (K(l)-K(8), KV, Al, A2) in COMl. The function 
returns a -2 if any format discrepancies are discovered, a 0 
if all parameters are correctly input. 

Usage: IERR = IPARAM (CDNO, JFMT, NO) 

CDNO - index to packed input card 
JFMT - Expected parameter format 
JFMT = 0 - Integer 
JFMT = 1 - Floating point 
NO. - Number of consecutive parameters 
to be checked 
IERR - Parameter check 

IERR = 0 - All parameters input correctly 
IERR - -2 - Format error discovered 


Common block COMl parameters: 


Common 

name 

Local 

name 

CARD (13) 

CAJ2P - Packed card image ^ '* : 

IN (5) 

IN - Input disk files 

L 

L - Number of lines 

C 

C - Number of columns 

K1-K8 

K (9) - Integer parameters 

KV 


Al 

A12 (2) - Real parameters 

A2 



Subroutine : UNPACK 

System library routines: CTOBIN 


Storage : 


365 


(8 


words 
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ISOPLT (Isodata Plot ) 

Purposes ISOPLT produces an isodata plot of the gray values from an 

input image. 

Usage: CALL ISOPLT (SMOOTH) 

SMOOTH - Option flag 

SMOOTH = .FALSE. - No spatial integration 
SMOOTH = .TRUE. - Spatial integration (averaging) 

is performed on each pixel of 
input image. 


Common block COMl parameters: 


Common 

Local 

name 

name 

IN (-5) 

I DP - Input file name 

ON 

IDO - Output file name 

L 

NL - Number of lines 

C 

NC - Number of colums 
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Method: 


i(7) 

i (8) 

i(l) 

i(6) 

X 

i (2) 

i (5) 

i (4) 

i (3) 


The above is an example of the numbering scheme used in the 
iso-data contouring routine for identifying the eight neighbors 
(i) of element x. 

The entire picture is smoothed by replacing each neighbor 
with the average of it and its eight neighbors, i.e., 

8 

x = x+ £ i (n) 

n-1 

9 

Each of the eight neighbors i (1 through "8) are searched to find 
at least one which is in the same contour plane or in a higher 
level plane as the center element (x) . 

Each of the four neighbors (2, 4, 6, and 8) are searched to 
find at least one which is in a lower contour plane than the 
center element (x) * 

Elements for which the neighbor conditions stated above are found to be true are 
marked as border elements and all others are set to zero. 

Subroutines: DATAIN, DISKR, DISKW 

Storage: 20774. words 


Required by: 


ISOGRAM 
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KONVL (Convolve) 

Purpose: KONVL applies a 3 by 3 moving window filter to an image 

by numerically convolving the input array with the 
3 by 3 filter. 

Usage: CALL KONVL 

Common block COM1 parameters: 


Common 

name 

Local 

name 

IN (5) 

JIN - Input file name 

ON 

JOUT - Output file name 


NLIN - Number of lines 

C 

NCLM - Number of columns 


Method: With filter W (J,K), input array I(C,L) and output array 

O (C, L) where C and L are the column and line coordinates 
respectively : 

3 3 

0 (M,N) = E E W(j,k) I(m-2+j, n-2+k) 




j-1 k-1 

Subroutines: 

DATAIN, 

DISKR, DISKW 

Storage : 

20407(8 

words 

Required by: 

filter 
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LABL (Label ) 

Purpose: LABL inserts labels, which are read from input cards, into 

an image. 

Usage: CALL LABL 


Common block C0M1 parameters: 


Common 

name 

Local 

name 

L 

NL - Number of lines in output frame 

c 

C - Number of columns (must be 1024) 

Kl 

NP - Number of label lines to be inserted 

K2 

JIN - input tape 

K3 

JOUT - Output tape 


Method: A packed display format tape is read down to the line at 

which the label is to be inserted. The label is overlayed 
and the process repeated until all labels have been inserted. 
The format of the label specification is: 

Col. Col. Col. Col. Col. 



J - Number of characters in the label 
K - Starting line position on the image 
L - Starting column position 
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The label begins in column 25 and continues for J columns. 

If the label is longer than 48 characters, it is continued 
on a second card beginning in column 1. 

Each character is a 9x7 matrix with a 9x3 matrix of trailing 
blanks; therefore, each character occupies 9 lines and 10 
columns. Line spacing is controlled by the user - a minimum 
of 5 lines between labels is recommended. 

The character set is: 

0 through 9, = ’ - () . , * / , A through Z 

The label is added to the original image so that the labels 
essentially overlay the image thus allowing overprinting, if 
desired. 

Subroutines : DSPLIN 

System library routines: KEDTPR, WRITER 

Storage: 3724 words 


Required by: 


LABEL 



15 September 


Purpose: 


Usage: 


1973 


69 System Development Corporation 

TM- (L) -HU-033/008/00 


MESS (Message) 


MESS is used by the diagnostics package to store and 
output specific error diagnostics. 

CALL MESS (FLAG) 

FLAG - Message indicator 


Storage: 


426 words 
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MLTAVS (Multi -Copy Average) 

Purpose: MLTAVS computes the average of up to five images. 

Usage: CALL MLTAVS 


Common block COM1 parameters: 


Common 

name 

Local 

name 

IN (5) 

T - Input file names 

ON 

0 - Output file name 

L 

L - Number of lines 

C 

M - Number of columns 

Kl 

S - Number of images to be averaged 


Method: For corresponding pixels of the images to be averaged , 

the mean gray scale value is calculated. The standard 
deviation is computed as 



Those gray scale values outside 1 standard deviation are elim- 
inated and a new mean value is calculated using only the gray 
scale values within 1 standard deviation of the original mean. 

Restrictions: (1) A maximum of 5 images may be averaged 

(2) The maximum size of the images to be averaged is 
2048 by 2048. 
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Subroutines: 


Storage : 


DISKR, DISKW 
3014 words 


Required by: 


AVERAGE 
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MPTORI (M/P to R/I) 

Purpose: MPTORI contains two entry points: 

MPTORI - Converts magnitude/phase complex file 
(interleaved line by line) to a real/ 
imaginary complex file (interleaved line 
by line) . 

RITOMP - Converts from a R/I file to a M/P file. 

Usage: CALL MPTORI 

CALL RITOMP 


Common block COM1 parameters: 


Common 

name 

Local 

name 


IN (5) 

IN - Input file name 

. 

ON 

10 - Output file name 


L 

NL - Number of lines 


C 

NC - Number of columns 



Restrictions: The number of columns (NC) must be less than or equal to 1024. 

- 4 1 

Subroutines: DISKR, DISKW 

System library routines: SIN, COS, SQRT, AT AN 2 

Storage: 10400. words 

to 


Required by: 


M/P TO R/I, R/I TO M/P 
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Purpose: 


Usage: 


MULDSP (Multiple Display Format) 


MULDSP arranges multiple images into a single frame in 
display format* Options include the expansion of each image 
by an expansion factor , and labelling under each row of 
images* 


CALL MULDSP (INTP, OUTP) 

IN TP - Input tape logical unit number 
OUTP - Output tape logical unit number 


Common block COM1 parameters : 


Common 

name 

Local 

name 

SN (5) 

SCRACH - Scratch disk file name 

L 

NL - Number of lines 

c 

NC - Number of columns 

K1 

NTP - Total number of images to be on output frame 


NTP <0 - no white spaces, butt images together 


NTP >0 - 100 lines of blanks between rows 

K2 

NPR - Number of images in one row across frame 

K3 

IOP - Option to invert gray scale values before output 


I OP - 1 - invert 


IOP - Other - No inversion 

K4 

I EXP - Expansion factor (a positive integer) 

K5 

NLAB - Total number of label cards following operator card 


Method: The images are taken from successive files on the input tape, 

expanded if required, and formatted into a single n x 1024 
frame in display format, where n is the variable length 
made up of NL * (Number of rows +■ 100) . Designated labels 
are written underneath the rows of images* A maximum of 
102 characters may be written per line. 
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Restrictions : 
Subroutines : 


System library 


Storage : 


All input images must be the same size. 
DISKR, DISKW, PL6 , FRAMWR, DSPLN 

routines: WRI TE R 

11161 words 

\o 


Required by: 


MULTIPLE DISPLAY 
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NOSFUG (Won-Symmetrical Function Generator) 

NOSFUG produces a non-symmetrical two dimensional function 
for use as a point spread function in deconvolution or for 
general filter operations where non-symmetry is necessary. 

CALL NOSFUG 


Common block COM1 parameters : 


Common 

name 

Local 

name 

ON 

JOUT - Output file name 

L 

N - Number of lines 

C 

M - Number of columns 


Purpose : 


Usage: 


Method: The user specifies the non-symmetrical function base in polar 

coordinates and the function profile as a function of the 
distance along any radial through the Z-axis extending to 
the base boundary. The base boundary is defined in polar 
coordinates (R and Theta). The function profile is defined 
as a percentage of any radial. 

■ 4 - 

Subroutines: DATAIN, DISKW 

System library routines: SQRT, AT AN 2 , LAGRAN 

Storage: 4360^g words 

Required by: NON-SYMMETRICAL PSF 
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NSA (Non-Symmetrical Alteration) 

Purpose: NSA performs a non-symmetrical alteration of the gray scale 

of an image based on correction planes which are input on 
tape unit A5 (8) . 

Usage: CALL NS A 


Common block COM1 parameters: 


Common 

name 

Local 

name 

IN (5) 

IN - Input file name 

ON 

ON - Output file name 

L 

NL - Number of lines 

C 

NC - Number of columns 

KV 

NP - Number of planes (3 - NP - 8) 


Input Data: The actual value of each correction plane is input following 

the operator card in table form as follows: 

CURVE=5 . , 10., 2(T. ,32., 45., 60.$ 

Tape Format: The values in CURVE represent the actual value of the gray 

schle of each plane. For example, if six planes are to be 
used for correction, they may be generated by photographing 
6 uniform fields of gray scales 5, 10, 20,... as indicated 
above. The digitized images resulting from photographing 
the uniform fields are the planes used for correction and 
are formatted on tape as follows: 
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Method : 


Restrictions : 


Subroutines : 
System library 


RECORD 1 
RECORD 2 


Line 1, Plane 1 
Line 1, Plane 2 




> 


RECORD NP 
RECORD NP+1 


Line 1, Plane NP 
Line 2, Plane 1 


• ^ 


RECORD NL*NP Line NL, Plane NP 


The tape, with the correction planes properly formatted, 
must be generated prior to operation. 


A three-point Lagrange interpolation is performed to find 
the corrected gray scale value for each pixel in the image. 


Assume the correction planes are designated X^, ^2* , *‘^n' 

The corrected value, GCX^) is then computed by interpolation 

as follows: Given the dependent variable table, X^^, 

X . . . .X , „ , as a function of the independent variable 

2(ij) n(ij) 

table, CURVE, find a value, G(X_), as a function of given 
argument, X^ . 


1) 3 - NP - 8 

2) Number of columns in image must be no greater than 1024. 

3) Number of columns in correction planes must be equal to 
number of columns in image. 


DISKR, DISKW, DATAIN 


routines : ' LAGRAN , EXIT 


Storage: 


22406 words 
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Required by: 


FIELD CORRECTION 
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ORAND (OR/ AND) 


Purpose: ORAND performs a logical AND or an inclusive OR on two input 

images , 

Usage: CALL ORAND (KO) 

KO - Option flag 

KO = 1 - Inclusive OR 
KO = 2 - AND 


Common block COM1 parameters: 


Coinmon 

name 

Local 

name 

IN(5> 

IN - Input file name 

ON 

ON - Output file name 

L 

NL - Number of lines 

C 

NC - Number of columns 


Subroutines : 

DISKR, DISKW, 

Storage: 

10222 , „ words 

Required by: 

AND', INCLUSIVE OR 
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Purpose: 

Usage: 

System library 

Storage: 
Required by: 


PACK (Pack 6 Characters Per Word) 


PACK transfers N characters which are stored as the least 
significant 6 bits of N words into an array where the 
characters are packed 6 per word. 

CALL PACK (A, N, B) 

A - Array of characters to be packed 
N - Number of characters 

B - Array into which characters will be packed 
routines : FLD 

130 /o words 

ALTER, DEPENDENT ALTER, PSF GENERATOR, FEATHER 
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PDGSA (Position Dependent Gray Scale Alteration) 

Purpose: PDGSA performs a position dependent alteration of the gray 

scale value of each pixel of an input image. 

Usage: CALL PDGSA 

Modification tables A and B, of up to 100 values each, 
represent the coefficients of the linear expression A + BX. 

Common block COMl parameters: 


Common 

name 

Local 

name 

IN (5) 

I - Input file name 

ON 

0 - Output file name 

L 

L - Number of lines 

C 

C - Number of columns 

KV 



N - Length of modification tables A and B 

i 


Method: The radial distance of a pixel from the center of the image 

is computed and is used as the index (m) to tables A and B in 
calculating a modified value as 

Modified value = A + B (original value) 

m m 

The modified value replaces the original value in the image. 

Restrictions: All alterations are symmetrical about the center point of 

the image. 
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Subroutines: DATAIN, DISKW, DISKR 

System library routines: SQRT 

Storage: 20342, words 


Required by: 


DEPENDENT ALTER 
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Purpose: 


Usage: 


Subroutines: 


PRERUN (Pre-Run Diagnostics) 

PRERUN checks an input card image for the presence of an 
equal - sign (to indicate a table) or a blank in card 
column 1 (to indicate a continuation card) . This subroutine 
is only called in the case that sense switch 3 is down 
(signifying a pre-run diagnostic run) and the scan routine 
in entry RDC of subroutine INIT has not found a legal operator 
name in columns 1-3. 

CALL PRERUN (CARD, 13 C) 

CARD - Packed (6 ch/wd) card image 
13C - Card columns 1-3 

UNPACK 


Storage: 


lOl^g words 
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PSFGEN (Symmetrical Point Spread Function Generator) 

Purpose: PSFGEN accepts a curve which defines a point spread function 

half-profile and creates a three-dimensional PSF array in 
floating point format. 

Usage: CALL PSFGEN 

The desired set,, of up to 50 floating point values 
representing the profile of the PSF to be generated, is 
input in tabular format . 


Common block COM1 parameters: 


Common 

Local 

name 

name 

ON 

JOUT - Output file name 

KV 

KSIZ - Size of the PSF to be generated 


Method : The three-dimensional PSF is generated by rotating the input 

curve about its left-hand axis. The distance from the center 
is calculated for each point and a gray value is calculated 
by interpolating between the input curve points using library 

routine LAGRNG. 

- J ' ’• 

Subroutines: DATAIN, DISKW 

System library routines: SQRT, LAGRNG 

Storage: 24166 ^g words 


Required by: 


PSF GENERATOR 
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RESTOR (Restore Disk File Directory) 

Purpose: RESTOR restores the copy of the file directory which was 

saved on disk by operator DUMPFD. 

Usage: CALL RESTOR (LOF, NF, DISK) 

LOF - Last opened output file 
NF - Operator parameter table 


Common block DISKFE 

DISK - Restore flag 

DISK = i TRUE . - Restore from disk 
DISK = .FALSE. - Restore from tape 

parameters : 

Common 

Local 


name 

name 


NLT (100) 

ITKO (464) 

' 

NCT (100) 



KST (100) 
NLPT (100) 


> Copy of Disk Track 0 

KDT (100) 

) 


KP (100) 
KUT (100) 

ITKl (464) 


IDT (100) 

\ 

> Copy of Disk Track 1 

NAF 



NOB 

LTU (2) 

NT (100) 

. ITKADR (7) 

J 



System library routines: DTAKE , DCHECK, DPUT 

Storage: 4211 (8 words 


Required by: 


RESTORE 
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REWND (Rewind) 

Purpose: REWND rewinds and optionally unloads a tape drive. 

The tape logical unit number and load/unload flag 
are parameters read from the IDAPS operator card. 

Usage: CALL REWND (CARD) 

CARD - Packed (6 ch/wd) representation of the 
operator card image 

System library routines: CTOBIN , CLOSE 

Storage: 66 worc * s 

Required by: REWIND 
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RSLGEN (Bar Chart Generator) 

Purpose: RSLGEN produces a two-dimensional bar chart 

consisting of alternating groups of black and white 
bars with spatial frequencies which range from one 
cycle/line to one-half cycle per column. 

Usage: 'CALL RSLGEN 

Common block COM1 parameters: 


Subroutines : 

DISKW 


Storage: 

4353 /£> words 
to 


Required by: 

BAR CHART 

""" ~ 








Local 

name 

JOUT - Output file name 
SIZE - size of bar chart 


Common 

name 

ON 

KV 
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SCAL (Scale ) 

Purpose: SCAL scales an input image to values between 0. and 63. 

Usage: CALL SCAL 


Common block C0M1 parameters: 


Common 

name 

Local 

name 

IN (5) 

IN - Input file name 

SN (5) 

IS - Scratch file name 

ON 

10 - Output file name 

L 

NL - Number of lines 

C 

NC - Number of columns 

K1 

KF - Curve flag 


KF = 1 - Logrithmic (base e) 


KF = 2 - Square root 


KF = 3 - Cube root 


KF = 4 - Square 


KF = Other - Linear 

K2 

LF - Clipping flag 


LF = 1 - Clip outside of BC and TC 


.LF = Other - No clipping 



A1 

BC — Lower clipping level. All values below BC will be set 


to zero. 

A2 

TC - Upper clipping level. All values above TC will be set 


to 63. 
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Method ; 


Subroutines : 


System library routines 


Storage: 
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The minimum (MIN) and maximum (MAX) value is determined 
by comparing each pixel in the image. 

The spread is computed as 
SPREAD = MAX - MIN 
and the step size, or increment, is 
STEP = SPREAD/64. 

The file is then rescanned to scale each pixel between 0. 
and 63. by the desired curve, 

DISKR, DISKW 
: ALOG , SQRT 

4476 words 


Required by: 


SCALE 
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Purpose: 


Usage: 


SCK (System Check) 

SCK is the main program in the disgnostic package. It 
contains individual checkout sections for each operator. 
General types of diagnostics made are: 

1) Check for mandatory non-zero, non-blank parameters 

2) Check size restrictions on disk files 

3) Check unique restrictions of each operator 

CALL SCK (ID, *, *,) 

ID - Operator index 

* - statement numbers in subroutine INIT, 

proceeded by $ , to which RETURN 1 and 
RETURN 2 refer in subroutine SCK. 


Common block COMl parameters: (See Appendix A) 

Subroutines: MESS, SZCK 


System library routines: . SQRT, SSWTCH, NOW 


Storage: 


2334, „ words 
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SEGMNT (Segment Subsets) 


Purpose: 


SEGMNT investigates an input image and identifies groups 
of pixels within a common subset. If a point within one 
subset is given, all other subsets are ignored. 


Usage: CALL SEGMNT 

Common block COMl parameters: 


Common 

name 

Local 

name 

i 

! 

i 

IN (5) \ 

JIN - 

i 

Input file names j 

SN (5) 

JSKT, 

JKST - Scratch file names . 

ON 

JOUT 

- Output file name ; 

L 

NL - 

Number of lines [ 

C 

NC - 

Number of columns \ 

Kl 

LL - 

Lower limit for thresholding . Pixels less than LL ; 

are excluded from subset calculation. t 

K2 

LU - 

Upper limit for thresholding. Pixels greater than 
LU are excluded from subset calculation. ; 

K3 

NULL 

- Size limit. Subsets smaller than NULL are not j 

included in the subset ID. j 

K4 

L^“ - 

Complement tag. j 
LN = 1 - Input file is complemented before processing.! 
LN = Other - Input file is not complemented. 

K5 

IL - 

Line coordinate of a point within a subset. 

K6 

i 

\ 

IC - 

Columns coordinate of a point within a subset . 


Method: The SEGMNT subroutine consists of five major sections: 

e Pixels which are members of a subset are determined by 

virtue of their being within the specified upper and lower 
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gray scale limits (LL and LU) and their having at least 
one neighbor within those limits* 

9 All members of a single subset are identified by assigning 

each pixel of the subset a unique identification number. 

If only one subset is desired (IL and IC > 0) , all other 
subset members are set to zero and the desired subset 
members are set to 63. 

<s> The number of subsets thus identified and the number of 

pixels which are members of each subset are counted and 
any subsets smaller than the specified size limit (NULL) 
are thrown out. 

© Subset numbers are reassigned to reflect the relative 

size of all remaining subsets. 

© An output array of subset ID numbers is prepared, and the 

operator statistics are printed. 

The basic algorithm involved is that which identifies 
members of a single subset. Assuming that a pixel yO 
has been identified as being a member of a subset and it 
has neighbors in the following locations: 


Y7 

y8 

yi 

ye 

y'O 

y2 

ys 

y4 

y3 


which may or may not be within a subset, the following 
flow diagram describes the algorithm in general terms* 
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Subroutines : 

DISKR, 

Storage: 

22234 ^ 

Required by: 

SUBSET 


DISKW 
, words 


C£ 



15 September 1973 


95 


System Development Corporation 

TM-(L) -HU-033/ 008/00 


Purpose: 


SHADE (Position Invariant Gray Scale Alteration) 


SHADE modifies the gray scale values of an image by using 
the original gray value as an index to a table of modified 
gray values and replacing the original value with the 
modified one. 


Usage: CALL SHADE 

The desired modification set of 64 values is input in 
tabular format. 


Common block Coml parameters: 


Common 

name 

Local 

name 

IN (5) 

JIN 

- Input file name 

ON 

JOUT 

- Output file name 

L 

L - 

Number of lines 

C 

c - 

Number of columns 


Subroutines: DAT AIN, DISKR, DISKVJ 

Storage : 465 ^ -words 


Required by: 


ALTER 
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Purpose: 
Usage : 


SHAVE (Frame Edge Rounding) 

SHAVE rounds off the edges of an input image, 

CALL SHAVE 

A table of up to 25 real values, representing percentages 
by which the elements on the edge of the frame will be 
modified, is entered in tabular format. 


Common block C0M1 parameters: 


Common 

name 

Local 

name 

IN (5) 

IN - Input file name 

ON 

10 - Output file name 

L 

NL - Number of lines 

C 

NC - Number of columns 

KV 

NP - Number of points in modification table 


Subroutines : 


datain, diskr, diskw 


Storage : 


$S£1 words 
to 


Required by: 


FEATHER 
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SHRINK 


Purpose : 


SHRINK reduces the size of an image by an integer factor. 


Usage: CALL SHRINK 

Common block C0M1 parameters: 


Common 

name 

Local 

name 

IN (5) 

IN - Input file names 

ON 

ON - Output file name 

L 

L - Number of lines in the picture 

C 

C - Number of columns in the picture 

Kl 

IFACT - Shrinkage factor 

K2 

KOP - Shrinkage method 

KOP = 1 - Every IFACT' th point 

KOP = Other - Average of consecutive IFACT x IFACT 
arrays 


Method : 


Restriction: 
Subroutines : 


An image that is ( IF ^ CT ' ) as lar 5® as the original image is 
formed by first calculating the average of consecutive 
IF*J«:T x IFACT subarrays from the original image, and then 
inserting the average value into consecutive pixels 
in a smaller array. 

The maximum column size is 1024. 


DISKR, DISKW 


Storage : 


10360 words 
(8 


Required by: 


SHRINK 
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SLICE 


SLICE plots the gray level values found on a line between 
two points of an image* 

CALL SLICE 


Common block COMl parameters: 


Common 

name 

Local 

name — — 

IN (5) 

IN - Input file name 

L 

L - Number of lines 

C 

c “ Number of columns 

Kl 

IY1 - Line coordinate point 1 

K2 

IXl - Column coordinate point 1 

K3 

IY2 - Line coordinate point 2 

K4 

1X2 - Column coordinate point 2 


Purpose : 


Usage: 


Method : 


Subroutines : 


The point slope form of the straight line equation is used 
to find the index of points found on a line between two 
points * 


DISKR 


Storage: 


10624 words 


Required by: 


SLICE 
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Purpose : 


Usage: 


SP1NT (Spatial Integration) 

SPINT performs a spatial integration of an image by 
averaging the eight immediate neighbors of each element 
and adjusting the element by some percentage. 


CALL SPINT 


Common block C0M1 parameters: 


Common 

name 

Local 

name 

IN (5) 

JIN - Input file name 

ON 

JOUT - Output file name 

L 

NL - Number of lines 

C 

NC - Number of columns 

Kl 

LMT - Allowable variance within which an element will 


not be adjusted 

K2 

KDT - Data type 

Al 

AJST - Correction factor (stated as a percentage) 


Method : 


' ''Eaph element E of the frame of data within a 1 element 
border is surrounded by 8 immediate neighbors N; i.e./ 


N 

N 

N 

N 

E 

N 

N 

N 

N 


The mean (X) of the 8 neighboring elements is computed for 
each element within the frame. If the center element is out- 
side a specified limit of the mean of its neighbors, it 
is replaced by adjusting it as follows: 
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Restriction: 


Subroutines: 


Storage : 


New E = E - [E - X) . AJST] 

For example, if E = 20, X= 28, LMT = 2, and AJST = ,8 {80%) : 

New E = 20 - [ (20 - 28) (.8)] 

= 20 - (-6.4) 

= 26.4 = 26 

The outer border of the image is unaltered. 

DISKR, DISKW 
20526 words 


Required by: 


SMOOTH 
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SPLT (Split Data) 


Purpose: SPLT has two entries: 

SPLT - Splits complex data which is interleaved line- 
by-line onto real (magnitude) components or 
imaginary (phase) components. 

MERG - Merges the separated complex data line-by-line 
into one file. 


Usage: CALL SPLT 

CALL MERG 

Common block COM1 parameters: 


Common 

name 

Local 

name 

IN (5) 

IN - Input file name 

ON 

IO - Output file name 

L 

NL - Number of lines 

C 

NC - Number of columns 

KV 

NP - Split flag 


NP = 1 - Split imaginary (or phase) 


NP = Other - Split real (or magnitude) 


1 N 


Restriction: 
Subroutines : 


Storage: 


The number of columns must be less than or equal to 1024. 
diskr, DISKW 
2234 ... words 


Required by: 


SPLIT, MERGE 



15 September 1973 


102 


System Development Corporation 
TM- (L) -HU-03 3/008/00 


Purpose ; 
Usage: 


Method : 


Restriction: 


Subroutines : 
Storage : 


SPREDl 

SPRED1 performs the multiplication of two complex vectors. 

CALL SPRED1 (N, IP, X, Y, GR, GI, MOL) 

N - Size of the arrays 
IP - Name of PSF file 

X - Real part of iamge transform array 
Y - Imaginary part of iamge transform array 
GR - Real part of PSF transform array 
GI - Imaginary part of PSF transform array 
MOL - Line control 

Let I = image = X + j Y 
P = PSF = GR + j GI 
0 R = real part of resulting image 
0 = imaginary part of resulting image, 
then (I) (P) = (X + jY) (GR + j GI) = 

(X) (GR) - (Y) (GI) + j [ (Y) (GR) + (X)(GI)] 

0 R = (X) (GR) - (Y) (GI) 

O x = (Y) (GR) + (X) (GI) 

The arrays must be the same size and less than or equal to 
1024 . 

DISKR 

14 2. _ words 


Required by: 


CONVOLVE 
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SPRED2 


SPRED2 performs the division of two complex vectors while 
allowing the limitation of the amount of magnitude change 
or phase change of the numerator resulting from this 
operation. 


CALL SPRED2 (NC, IN, IMR, IMI, PSFR, PSFI, NMLIM, NPLIM, MOL ) 

NC - Number of columns 
IN - PSF file name 

IMR - Real part of the numerator (image transform) 
array 

IMI - Imaginary part of the numerator array 
PSFR - Real part of the denominator (PSF transform) 
array 

PSFI - Imaginary part of denominator array 
NMLIM - Limit on magnitude amplificaton (integer) 
NMLIM = Blank - No limit 

NMLIM =0 - Allow no magnitude change 

NMLIM = Other - Limit all amplifications 

greater than NMLIM to NMLIM 
NPLIM - Limit on phase amplification (integer) 

NPLIM = Blank - No limit 

■\^ ' NPLIM = 0 - Allow no phase change 

NPLIM = other - Limit all phase changes 

whose magnitude is greater 
than NPLIM to NPLIM. 

MOL - Line control 

The real/imaginary components are converted into magnitude/ 
phase components. 

Let A = image magnitude vector 
0 = image phase vector 
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B = PSF magnitude vector 

0 = PSF phase vector 

C = resulting magnitude vector 

a = resulting magnitude vector 

= A / 8-0 
B 


. A/e 
c /a = - — 

B/JET 


If NMLIM 


If NMLIN 
If NMLIN 


If NPLIM 

If NPLIM 
If NPLIM 


Blank - operation above is carried out 

for all PSF points not equal to 

zero. If a PSF point is equal 

zero, the corresponding point of 

the image magnitude vector is unaltered. 

0 - each point of B is set to one. 

Other - for those points such that 

1/B^ is greater than NMLIN, set 

1/B . to NMLIN. 
l 

Blank - operation above is carried out 
completely. 

0 - 0 is set to zero. 

Other and | 0 | is greater than MPLIM , 

0 is set to sign ( 0 ) times NPLIM. 


Restriction: 


Subroutines: 


The arrays must be the same size and less than or equal to 
1024. 


DISKR 


System library routines: CTOBIN, SQRT , ATAN2, COS, SIN 


Storage: 


450 (8 


words 


Required by: 


DECONVOLVE 
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Purpose: 

Usage: 


Method : 


Restrictions : 
Subroutines : 
Storage: 


SPRED3 


SPRED3 performs the multiplication of two complex vectors, 
and computes the complex conjugate of the result. 

CALL SPRED3 (N,IP,X,Y,GR,GI,M0L) 

N - Size of the arrays 
IP - Name of PSF file 

X - Real part of image transform- array 
Y - Imaginary part of image transform array 
GR - Real part of PSF transform array 
GI - Imaginary part of PSF transform array 
MOL - Line control 

Let I = image = X + j Y 
P = PSF = GR + j GI 
0 *= real part of resulting image 

Oj. = imaginary part of resulting image 
then (I) (P) = (X + jY) (GR + j GI) = 

(X) (GR) - (Y) (GI) + j [ (Y) (GR) + (X) (GI)] 

0 R = (X) (GR) - (Y) (GI) 

0 I = -(Y) (GR) - (X) (GI) 

The arrays must be the same size and less than or equal to 1024. 
DISKR * 

142, 0 words 


Required by: 


CORRELATE 
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STAT (Statistics ) 

Purpose: STAT produces a histogram of the distribution of the gray 

scale values of an image and computes the following statistics: 

Sum of gray scale values 
Mean gray scale value 
Average deviation 
Variance 

Standard deviation 

Isolated maximum and minimum points 
Usage: CALL STAT 


Common Block COMl Parameters : 



Subroutines: DISKR 

^'V k - * 

System library routines: ^ SQRT 


Storage: 5063, „ words 

1.0 


Required by: 


HISTOGRAM 
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Purpose: 

Usage: 


Note: 


STKPIC (Stack Frames of Data on Tape) 

STKPIC transfers an image from disk and stacks it on tape. 

CALL STKPIC (10) 

10 - Output tape logical unit number 

The output tape, 10, is not rewound, and no end of file is written 
after the last record. 


Common block COM1 parameters: 


Common 

Local 

name 

name 

IN 

IN - Input file name 

L 

NL - Number of lines 

c ; 

NC - Number of columns 


Subroutines: 


Storage: 
EEquired by: 


DISKR 

1153 words 
STACK 
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Purpose : 


Usage: 


SZCK (Size Check ) 

SZCK is used by the diagnostic package to provide size 
information about existing input files (found in the disk 
file directory). 

CALL SZCK (IN, ISIZE) 

IN - Input file array 

ISIZE - Array containing the number of lines and 
number of columns of each file in IN. 


Storage: 


1716„. words 
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TFGEN (Transfer Function Generator) 

Purpose: TFGEN generates a transfer function from user-provided 

independent variables. 

Usage: CALL TFGEN (A,B ,C ,D) 

A - Independent variable 
B - Independent variable 
C - Independent variable 
D - Independent variable 

Common block COMl parameters : 


Common 

name 

Local 

name 

CARD (13) 

CARD - Input card image (packed) 

ON 

ON - Output disk file name 

L 

NL - Number of lines 

C 

NC - Number of columns 

K1 

K1 - Transfer function code 


Kl = 1 - Polynomial 


Kl = 2 - Delta function 


Kl = 3 - Gaussian 


^ Kl - 4 - MTF of defocussed perfect lens 


Kl = 5 - Quarter wave of zero order Bessel function 


Method: Based on the basic transfer function profile specified by the 

user and any required independent variables that must be 
provided for the selected profile, a one dimensional table is 
generated from the appropriate mathematical expression. The 
table length is selected equal to the number of picture 
elements from the center of the two dimensional output array 
to its outer edge. 
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The function profile curve is then "rotated" about the 
center point to generate a two dimensional output array 
having such rotational symmetry that any plane passing 
through the center point, normal to the base, will contain 
the original one dimensional profile curve and its mirror 
image. This is done by calculating the distance of every 
picture element in the output array from the center element 
and then finding the value for that element on the profile 
curve. A Langrangian interpolation technique is used to find 
values for those point-to-center distances which are not an 
integer number of picture elements* 

The two dimensional array of real data values is written as an 
output file and the values of the one dimensional profile from 
which the array was produced are printed on the on-line printer. 


The transfer functions T(x) have the following mathematical 
formulae — 

K = 1 T(x) = A + Bx + Cx 2 + Dx^ 

K = 2 T(x) = A where line = B 

and column = C 

= o where line ^ B 

and column / C 

(If no B and C are specified, the center point = A) 

X l/3x. 2 

„ „ 2 'A ’ 

K - 3 T(x) = Be 



K = 5 
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Subroutines : 
System library 
Storage: 
Required by: 


DISKW, BESJ 

routines: SQRT, ATAN2 , LAGRAN 

11046 words 
MTF 
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Purpose : 

Usage: 

System library 
Storage : 


TIME 


TIME prints out the time required to execute the preceding 
operator, if sense switch 5 is not down. 


CALL TIME (KSS) 

KSS - Time at beginning of operator's execution. 


routines: SSWTCH, NOW 


72 


(8 


words 
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Purpose : 


Usage: 


Subroutines : 
System library 
Storage: 
Required by: 


TRANS (Transpose) 

TRANS performs a transposition, and associated operations^ 
on a core buffer, and writes the final image out on the 
output file (A6) . 


CALL TRANS (OUTB, COREB, OUTF, NL, NUM, F, MAXC, IT, INDEX, 
REVOLV) 

OUTB - Output buffer 

COREB - Core buffer 

OUTF - Output file on disk 

NL - Number of lines in the image 

NUM - Number of columns to transpose as rows 

F - Format 

F = 1 - Integer 
F = 0 - Floating Point 

MAXC - Maximum number of lines to write (actual number 
of unpacked characters in original line) 

IT - Current line on output disk file. If negative, 
then a horizontal rotation has been selected. 
INDEX - Variable dimension indicator for COREB 
REVOLV - Transposition option indicator 

REVOLV = 1 - Transpose only. (IT = 1). 

REVOLV = 2 - Transpose and vertical rotate. (IT 
„ REVOLV = 3 - Transpose and horizontal rotate. 

(IT = -NC). 

REVOLV = 4 - Transpose and both. (IT = -NC) . 


= 1 ). 


FLIP, DISKW 

routines : FLD 

363 , Q words 
(o 

TRANSPOSE, +90 DEGREE ROTATE, -90 DEGREE ROTATE, MINOR 
TRANSPOSE, VERTICAL ROTATE, HORIZONTAL ROTATE, 180 DEGREE ROTATE 
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TRNP5 (Transposition Control) 

TRNPS controls the rotations required by the following 
operators : 

1) TRANSPOSE 

2) +90 DEGREE ROTATE 

3) -90 DEGREE ROTATE 

4) MINOR TRANSPOSE 

5) VERTICAL ROTATE 

6) HORIZONTAL ROTATE 

7) 180 DEGREE ROTATE 

CALL TRNPS (COREFL, INDEX, REVOLV, NL) 

COREFL - Buffer of core storage 

INDEX - Variable which dynamically dimensions COREFL 
to maximum advantage 
REVOLV - Rotation options Indicator 

REVOLV = 1 - Transpose (used for TRANSPOSE 
operator) 

REVOLV = 2 - Transpose and Vertical Rotate 

(used for +90 DEGREE ROTATE operator) 
REVOLV - 3 - Transpose and Horizontal Rotate 

(used for -90 DEGREE ROTATE operator) 
v ' >. REVOLV = 4 - Transpose and both Vertical and 

' ■ Horizontal Rotate (used for MINOR 

TRANSPOSE operator) 

REVOLV = 5 - VERTICAL ROTATE 
REVOLV = 6 - HORIZONTAL ROTATE 

REV0LF = 7 - Both Vertical and Horizontal Rotate 

(used for 180 DEGREE ROTATE operator) 
NL - Number of lines in input image 
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Common block C0M1 parameters: 


1 Common 
name 

Local 

name 

IN(5) 

INFILE - Input file name 

SN(5) 

SCRACH - Scratch file names 

ON 

0UTFIL - Output file names 

L 

NL - Number of lines 

C 

NC - Number of columns 

K1 

FMT - Format of the data 
FMT = 1 - Integer 
FMT = Other - Floating point 


Method: The input image is stored in COREFL and up to 5 scratch files. 

The core memory portion is transposed and written on the output 
file. Each scratch file is transferred to memory, one core 
buffer section at a time, and transposed and output in the 
same way. If possible, the scratch files are dimensioned the 
same size, or a multiple of the size, of the core buffer to 
minimize disk transfer time. 


Restrictions: The maximum image size is 2048 x 2048. 

Subroutines: DISKR, DISKW, FLIP, PL6 , TRANS 


Storage: 3145 words 

Required by: TRANSPOSE, +90 DEGREE ROTATE, -90 DEGREE ROTATE, MINOR 

TRANSPOSE, VERTICAL ROTATE, HORIZONTAL ROTATE, 180 DEGREE 
ROTATE 
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Purpose : 
Usage: 


Method : 


Storage: 

Required 


TSFSZ (Transpose Scratch File Size) 

TSFSZ calculates the optimum sizes for the scratch files used 
by the TRANSPOSE and other related operators. 

CALL TSFSZ (NL, NC, FMT, NO, KFN) 

NL - Number of lines in original image 
NC - Number of columns in original image 
FMT - Format of original image 
FMT = 1 - Integer 
FMT = 0 - Floating point 

NO - Dummy integer. Returns number of scratch files to 
be opened. 

KFN - Dummy array. Returns size (number of columns) 
calculated for each scratch file. 

The optimum size of the scratch files to be opened is deter- 
mined by the following considerations: 

1) There may be a maximum of 5 scratch files opened. 

2) The disk files should fit into the core buffer, with 
one core block per disk file, if possible. If not, 
the size of the disk file should be a multiple of the 
size of the core buffer (10240). 

3) If unequal sizes are necessary between the five files, 
the last, or fifth, file should contain the overflow. 

272 , c words 
(.o 

by: TRANSPOSE, +90 DEGREE ROTATE, -90 DEGREE ROTATE, MINOR 

TRANSPOSE, VERTICAL ROTATE, HORIZONTAL ROTATE, 180 DEGREE 
ROTATE 
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Purpose : 

Usage: 

System library 
Storage: 


UNPAC (Unpack 1 character from packed array) 

UNPAC is a function which extracts the j th character from 
a packed array. 

B - UNPAC (A, J) 

A - Array of packed data 
J - Index of character to be unpacked 

B - Desired character right justified with leading zeros 
routines : FLD 

103 words 
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Purpose : 

Usage: 

System library 
Storage: 


UNPACK (Unpack 1 character from packed array) 

UNPACK is a function which extracts the character from 
a packed array. 

B = UNPACK (A , J ) 

A - Array of packed data 
J - Index of character to be unpacked 
B - Desired character, right justified with leading 
blanks 


routines: FLD 


103 (8 


words 


Required by: 


ALTER, DEPENDENT ALTER, PSF GENERATOR, FEATHER 
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Purpose : 


UPL6 (Unpack Line - 6 Characters/Word) 


UPL6 contains two entry points: 


UPL6 - Unpacks 6 bit characters (6 char/word to 1 char/ 
word) . 

PL6 - Packs 6 bit characters (1 char/word to 6 char/word). 


Usage: CALL UPL6 (LP, SC, NC, LU) 

CALL PL6 (LU, SC, NC , LP) 

LP - Packed line 

SC - Starting character to be packed /unpacked 
NC - Number of characters to be packed /unpacked 
LU - Unpacked line 

System library routines: FLD, PCKINT , UNPINT 

Storage: 543 words 
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XPAND (Frame Expansion) 

Purpose: XPAND performs an expansion or an extraction of a selected 

portion of an image. 

Usage: CALL XPAND (NP) 

NP - Expansion factor 

Common block C0M1 parameters : 


Common 

name 

Local 

name 

IN (5) 

■ JIN - Input file name 

ON 

JOUT - Output file name 

L 

L - Number of lines in original frame 

C 

C - Number of columns in original frame 

Kl 

LUL - Left upper line element number 

K2 

LUC - Left upper column element number 

K3 

NE - Number of lines in extracted portion of image 

K4 

NC - Number of columns in extracted portion of image 


Method: Each line of the original image is read from the input tape 

and v each element within the area to be expanded is duplicated 
NP times. Each line of data is then output NP times on the 
output tape. To extract rather than expand a selected 
portion of an image, NP should be set to 1. 

Restrictions: The product of the size of the area to be expanded or extracted 

and NP may not exceed 2048. 


Subroutines : 


DISKR, DISKW 
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System library routines: EXIT 

Storage: 10302, Q words 


Required by: 


EXPAND* EXTRACT 
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APPENDIX A 
COMMON BLOCKS 


COMMON/COMl/CARD,IN,SN,ON,PN 1 L,C,Kl,K2 ) K3,K4,K5,K6 t K7,K8,KV,Al,A2 

Common block C0M1 is used by the IDAPS system to store information used throughout 
the system. Although called by various local names, the information can be 
generally defined as follows: 


CARD (13) 
IN (5) 
SN<5) 

ON 

PN(5) 

L 

C 

K1-K8 

KV 

A1,A2 


- Input operator card image, packed six characters per word 

- Array of assigned input file names 

- Array of assigned scratch file names 

- Output file name 

- Temporary file names used for format conversion of input 
files 

- Number of lines, as defined by the last processed SIZE 
operator 

- Number of columns, as defined by the last processed SIZE 
operator 

- Integer parameters from operator card 
Integer parameter from operator card 

- Real parameters from operator card 


C0MM0N/DISKFD/NLT , NCT , KST , NLPT , KPT , KP ,KUT , IDT ,NAF ,N0B , LTU , NT , ITKADR 

V 

Common block DISKFD is used by the IDAPS system as a disk file directory 
containing information about all disk files (up to 100) which are opened. The 
directory is divided into arrays which can be defined as follows: 


NLT(100) 
NCT (100) 
KST (100) 
NLPT(IOO) 
KDT(100) 


- Number of lines in the image 

- Number of columns in the image 

- Starting track number 

- Number of lines per track 

- Data type 

KDT = 0 - Floating point 
KDT = 1 - 6-bit integer packed 
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KP(100) 


KUT(IOO) 
IDT (100) 
NAF 
NOB 
LTU (2) 

NT 

ITKADR(7) 


- Pointer to buffer 

KP = 0 - Closed 

KP = 1-7 - Buffer 1-7 

- Unit specification (0 or 1) 

- File identification 

- Number of assigned files 

- Number of open buffers 

- Last track used 

- Number of tracks 

- Track address of data currently in buffer KP 


COMMON/BREAK/KBREAK 

Common block BREAK is used to relay the information that the break key (sense 
switch 2) has been pressed by the user. 

COMMON / PRD / 1 ERR 

Common block PRD is used to relay the information that an error condition has 
been discovered by the diagnostic package. 


COMMON/TIMS /KSS 

Common block TIMS is used to store the starting time for an operator's execution, 
when the timing option (sense switch 5) is discovered. 

C0MM0N/C0M2/Y ,GR ,GI 

Common block COM2 is used to store intermediate results and to share common 
buffers between the Fourier filtering subroutines (FRX2V , CONV, and FFIL) . The 
names and dimensions defined for COM2 in the subroutines are as follows: 


FRX2V 

CONV 

FFIL 

Y(512) 

DUM(1536) 

FIPIC (512) 

GR(512) 


EIPIC (512) 

GI (512) 


FILTER (512) 
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APPENDIX B 
DIAGNOSTICS 


IDAPS 1 diagnostics package produces error messages which are output on the 
on-line printer. One group of these messages, primarily concerned with syntax 
errors, originates from the pre-run diagnostics subroutine, SCK, and is stored 
in subroutine MESS. Following printout of one of this group of messages, two 
paths of operation may be followed. If the system is in the pre-run diagnostic 
mode (sense switch #3 down) , the checkout of the remaining input cards is 
continued. During normal execution, however, the system pauses and awaits 
input of a corrected card before continuing execution. No execution of the 
incorrect operator is performed, any opened input disk files are closed, and 
output and scratch files for that operator are deleted. The following chart 
lists the non-fatal error messages in this group, the IDAPS operators to which 
they apply, and a brief explanation of each. 
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DIAGNOSTICS 


NON-BLANK PARAMETER REQUIRED 

TAPE LOGICAL UNIT NUMBER MUST BE 8, 10, or 11 

INPUT FILE PARAMETER INCORRECT 
WRONG NO. INPUT FILE NAMES SPECIFIED 


INPUT FILES MUST BE THE SAME SIZE 
INPUT FILES MUST BE SQUARE 

DIMENSION MUST BE A POWER OF 2 

DECONVOLVED IMAGE LESS THAN INPUT IMAGE 
PSF GREATER THAN INPUT IMAGE 


OPERATORS 


EXPLANATION 


A blank occurs where a number must be. 

LABEL Wrong number given for tape parameter. 

DUMP FILE 
LOAD FILE 

AVERAGE Number of files parameter is either too many 

LOAD FILE or missing. 

AVERAGE Number of files listed not equal to number of 

CONVOLVE files designated by a parameter, or by 

CORRELATE definition of the operator. 

MERGE 

MATH 

INCLUSIVE OR 

AND 

INSERT 

MOMENT j 

AVERAGE Lines and columns of two or more input files 

MOMENT are not equal. 

FOURIER FILTER 

DECONVOLVE Lines are not equal to columns in input file(s). 

CONVOLVE 

CORRELATE 

FOURIER FILTER 

DECONVOLVE Dimension of input file is not a power of 2 

BAR CHART between 32 and 512. 

CONVOLVE 

CORRELATE 

DECONVOLVE Parameter designated size of square decon- 

volved image is less than the number of lines 
in the input file. 

DECONVOLVE Parameter designated the size of the square 

PSF is greater than the number of lines in 
the input file. 
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DIAGNOSTICS 

LINES AND COL * OF 2ND FILE MUST = COL* OF 1ST 

N CANT BE LESS THAN HALF THE PICTURE DIAGONAL 
MOD CONTAINS MORE THAN 25 VALUES 

FILE MUST HAVE 1024 COLUMNS PER LINE 

NON-ZERO PARAMETER REQUIRED 
OPERATION FLAG OUT OF BOUNDS 


FRAMES ON ROW GREATER THAN TOTAL NO. FRAMES 

M> 

CM 

^ PSF IS GREATER THAN 100 


RESOLUTION CHART IS LARGER THAN 2048 
PARAMETER OUT OF BOUNDS 


LARGER FILE MUST PRECEDE SMALLER ONE 

SUBFRAME MUST FIT +2 WITHIN LARGER ARRAY 
SUBFRAME OVERLAPS EDGES OF LARGE ARRAY 

NO. LINES OF 1ST FILE MUST BE 2*NO* COLS* 

NO* COLUMNS MUST BE NO GREATER THAN 1024 


OPERATORS 


EXPLANATION 


FOURIER FILTER Number of lines of the square second input 

file does not equal the number of columns 
of the first input file* 

/ 2 2 

/ NL NC 

DEPENDENT ALTER Parameter N is less than / (— ) + (— ) 

FEATHER Parameter designated the number of values 

in table MOD is greater than 25. 

LABEL Input file has less than or more than 1024 

columns . 

SIZE A zero number of lines or columns is illegal. 

MATH Parameter is either less than 1 or greater 

than 8. 

MULTIPLE DISPLAY Parameter indicating the number of pictures/row 

must be less than the total number of pictures. 

PSF Parameter indicating the size of the PSF is 

greater than 100 or less than zero. 

BAR CHART Parameter indicating the size of the bar chart 

is greater than 2048* 

FIELD CORRECTION Parameter less than 3 or greater than 8 for 

SCALE FIELD CORRECTION. Minimum clipping value 

greater than maximum clipping value for SCALE. 

INSERT Second input file is larger than first input file. 

INSERT Size of second input file too big to allow inser- 

tion into first input file* 

FOURIER FILTER Incorrect size for input file. 

FIELD CORRECTION Number of columns of input file greater than 1024. 
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Special disk file diagnostic messages handle errors which occur when the disk 
is being read or written on and during the file directory bookkeeping. These 
messages generally indicate irrecoverable error conditions and an exit from 
the IDAPS system occurs. 

Other execution time diagnostics are not fatal. Either execution is continued 
following output of the error message, or instructions to the user are printed 
for re-entrance of the erroneous card. 

A list of these diagnostics, whether they are fatal or not, the subroutine 
from which they originate, and a brief explanation of each is given in the 
following chart. 
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SUBROUTINE EXPLANATION 

ASSIGN Operator xxx requires more than one input file 

and the names have not been provided in 
parentheses on the operator card. 

ASSIGN No comma or close parentheses was found for 

multiple input file names specified on 
operator card. 

ASSIGN A file name has occurred more than once in a 

multiple input file name specification. 

ASSIGN Either too many or too few names have been 

provided in the multiple input file name 
specification in parentheses. 

' Maximum image size has been exceeded for an 
integer array. 

An attempt has been made to open more than the 
maximum of 7 image files at one time. 

DISKD An attempt has been made to delete a non- 

existant file. 

DISKC An attempt has been made to close a file that 

has been previously deleted or one that has 
not been opened. 

DISKR An attempt has been made to read/write on a 

DISKW file which has not been opened (i.e., a buffer 

has not been assigned to it) . 

DISKR An attempt has been made to read/write on a 

DISKW file which does not have an entry in the disk 

file directory. 


DISKO 

DISKO 
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DIAGNOSTICS 

FATAL SUBROUTINE 

EXPLANATION 

ERROR - SUBSET NUMBER ASSIGNMENT HAS EXCEEDED 
AVAILABLE BUFFER SPACE 

SEGMNT 

More than 2048 subsets have been found which 
fall within the specified thresholds. 

TOO MANY CHARACTERS IN DISPLAY LINE - EXCESSIVE 
CHARACTERS IGNORED 

DSPLIN 

More than 102 characters provided for display. 

DATA CARD FORMAT ERROR - CLEAR CARD READER, 
CORRECT AND RELOAD CARDS, READY CARD READER 

DATAIN 

No equals sign was found on first input data 
table card. 

INCORRECT NUMBER OF DATA VALUES PROVIDED - 
CLEAR CARD READER, CORRECT AND RELOAD DATA 

DATAIN 

The count of the input values in the tables 
is not equal to the number of values expected. 

ONE OR MORE PARAMETERS IN WRONG FORMAT 

INIT 

No decimal point found where real parameter 
required, or decimal point found where integer 
parameter required. 

xxx IS AN ILLEGAL OPERATOR 

PRERUN 

, A card input during a pre-run diagnostic 
check is not part of a table or a continuation 
card and the first three characters are not 
a legal IDAPS operator. 

XXX NOT LEGAL OPERATOR - CLEAR CARD READER, 
CORRECT AND RELOAD CARDS, READY CARD READER 

RDC(INIT) 

First three characters of name on operator 
card (card columns 1-3) are not recognized 
as an IDAPS operator. 

ERROR IN INTERPOLATION ROUTINE 

FUNGN 

System error has occurred in system library 
routine 'LAGRNG. 

INTERPOLATION ERROR 

NSA 

System error occurred in system library inter- 
polation routine LAGRAN. 

ERROR - xxxxxxxx DATA TRANSFER MAY BE INCOM- 
PLETE OR WRONG 

DISKW 

System error xxxxxxxx has occurred while 
attempting to write to disk (system library 
routine DCHEK) . 

DATA TRANSFER ERROR — xxxxxxxx 

DUMPFD 

RESTOR 

System error xxxxxxxx has occurred while 
attempting to dump a disk file or write on 
disk (system library routine DCHEK) . 
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APPENDIX C 

PROCESS FOR ADDING AN OPERATOR TO THE IDAPS SYSTEM 
A. Subroutine IDAPS2 

The array NF must be altered to include information about the new operator. 
Add 

,6HXXX ,A,B,C/ 

to the DATA declaration for array NF, where XXX are the first three 
characters in the operator’s name, and A,B, and C are defined as follows: 

A — Number of input disk files to be assigned 
B - Number of scratch disk files to be assigned 
C - Number of output disk files to be assigned 

The DATA declaration for integer NO and the implied DO delimiter in the 
DATA statement for array NF , must be incremented to reflect the additional 
operator. 

The section of subroutine IDAPS2 pertaining to the new operator should 
be in the following format: 

nn CALL OPnn 

CALL SCK (I, $nn .00, $6 00) 
nnOO CONTINUE 

CALL NEW0P 
GO TO 600 

V*. 

where nn by convention indicates the new operators index into array 
NF (also contained in integer I) . OPnn is the entry for the operator in 
subroutine INIT, SCK is the driver for the diagnostic package, and 
NEW0P is the new execution subroutine or subroutine package* 

At statement number 200, the GOTO switch must contain the statement 
number nn. 
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B. Subroutine INIT 

The new operator, if it requires the opening of any disk files or the 
input of any parameters, must have an entry in subroutine INIT. Here, 
the parameters are read in by calls to IPARAM, format conversion is per- 
formed by a call on CKFMT, and any needed disk files are opened by calls 
to DISKO. 

C. Subroutine SCK 

Any diagnostics to be made on the new operator should be checked in 
subroutine SCK. As in IDAPS2, the operator index is used as a statement 
number and must be included in the initial GO TO switch. 

Checks on dimensions of disk files may be made by calls on SZCK. The 
diagnostic messages are stored in array B in subroutine MESS. Any 
additions made to this array must be reflected in the dimensions of the 
array. The maximum number of characters in a diagnostic message is 46. 

D. Subroutine "NEWOP" 

After its proper linkage, N, is determined, the subroutine (or subroutine 
package) NEWOP is inserted in the run deck in the following manner: 




